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Abstract  o£  Dissertation  Presented  to  the  Graduate  Council  In 
Partial  Fulfillment  o£  the  Requirements  for  the 
Degree  of  Doctor  of  Philosophy 


THE  EFFECTS  OF  MAJOR  IMPURITIES  ON  THE  TRANSPORT 
OF  HYmOGEN  THROUGH  PALLADIUM 


By 

Leon  Lee  Chiu 

April,  1966 

Chairman:   Prof.  Robert  D.  Walker,  Jr. 
Major  Department:   Chemical  Engineering 

The  purpose  of  the  investigation  was  to  study  the  effect  of 
major  impurities  on  the  transport  of  hydrogen  through  solid  palladium- 
silver  electrodes. 

Experimental  data  were  to  be  gathered  on  the  effect  of 
impurities  such  as  CO-,  CO,  N-  and  CH.  on  the  transport  rate  with 
temperature,  pressure  and  bleed  rate  as  parameters.  These  impurities 
might  be  expected  to  be  present  in  the  effluents  from  ammonia  or 
hydrocarbon  reformers.  The  temperature  range  covered  extended  from 
300°C  to  460°C  and  pressure  range  from  1069  mm  Hg  to  2310  mm  Hg. 

The  transport  rate  of  pure  hydrogen  through  0.004  inches  thick 
palladium-silver  foil  can  be  described  by  an  equation  of  the  form: 

1 

\   -  \^   exp  (^-  -  J  =  3.8  X  10  P    exp  [-  --^  J 

2 

where   N.  "   transport  rate  in  gm. moles/cm  sec. 


xiv 


t 


P  =»  pressure  differential  in  mm  Hg 

R  ■  gas  constant 

T  «»  absolute  temperature  in  **K 

AE  =  energy  of  activation 

«, ,~  ="  constant 
1  n 


Calculation  of  diffusivity  in  the  palladium- silver  alloy  gives 
an  equation  of  the  form: 

D  =  5.9  «  10-^xp  (-  ^  ) 

Comparison  of  this  diffusivity  with  that  for  pure  palladium  indicates 
that  introduction  of  silver  lowers  the  diffusivity. 

Introduction  of  impurities  lowers  the  transport  rate  of 
#^  hydrogen  through  palladium-silver  foil.  This  can  be  attributed  to 

the  following  effects:   (1)  hydrogen  is  diluted  by  the  impurities; 

(2)  impurities  hinder  the  diffusion  of  hydrogen  in  the  gas  phase; 

(3)  occupation  of  adsorption  and  dissociation  sites  by  impurities.  The 
overall  effect  of  impurities  is,  therefore,  to  lower  the  concentration 
of  hydrogen  atoms  in  the  upstream  surface  of  the  palladium-silver  foil, 
and  since  the  diffusion  rate  in  the  alloy  is  directly  proportional  to 
concentration,  the  transport  rate  decreases. 

^  The  transport  rate  of  hydrogen  through  0.004  inches  thick  25% 

palladium- silver  foil  with  impurities  can  be  described  by  the  following 
relationship  at  certain  temperature  and  pressure: 


"*  ■  -A„.x  ir^  =  "i^'e-H.  (-  i  ) 


B 


B-k2 


"^"^   *'Amax  "  °«''i'»"™  transport  rate  obtainable,  i.e.,  when  the  bleed 


rate  is  infinite 
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B  -  bleed  rate 

k.,k.  =>  constant 

1  .  , 

—  a  exponential  constant 

n 

^  AE  =  energy  of  activation 

Comparison  of  the  equation  for  N^^^  with  N  for  pure  hydrogen  shows 
that  the  value  of  --  is  lower  and  the  energy  of  activation  AE  is 
higher  when  impurities  are  present.  This  suggests  that  surface 
effects  are  present  when  impurities  are  introduced.  Furthermore,  the 
term  B/(B  -  k^)   is  probably  caused  by  effect  (2).   It  is  also  expected 
that  because  of  effect  (1),  the  partial  pressure  of  hydrogen  should 
be  used  in  the  equation  instead  of  total  pressure. 
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CHAPTER  I 
INTRODUCTION 
A.  General 


Since  Thomas  Graham  (1)  investigated  the  transport  o£  hydrogen 
through  palladium  in  1868,  numerous  investigators  (2,3,4)  have  made 
the  same  study.  However,  the  results  vary  from  one  investigator  to 
another,  sometime  as  much  as  a  thousand  times.  Recently,  some 
Investigators  (5,6,8)  have  succeeded  in  narrowing  this  variability  by 
simply  using  a  continuous  system.  Previous  investigators  almost 
exclusively  used  static  systems,  which  when  the  £eed  gas  contains 
small  amounts  of  impurities,  permit  the  impurities  to ;< accumulate  on;:the 
upstream  side  of  the  palladium  membrane,  thereby  producing  erroneous 
results.  A  continuous  system  bleeds  part  of  its  feed  gas  to  the 
atmosphere,  so  that  the  impurities  do  not  accumulate  on  the  upstream 
side  of  the  membrane. 

The  use  of  a  continuous  system  does  not  entirely  eliminate  the 
variability  of  the  results,  for  there  are  also  other  factors  involved. 
In  view  of  the  fact  that  the  transport  of  pure  hydrogen  through 
palladium  is  still  subject  to  controversy,  the  transport  of  pure 
^  hydrogen  has  been  subjected  to  detailed  analysis.  These  factors 

which  will  be  studied  are  the  thermodynamics  of  the  palladium- hydrogen 
system,  which  is  important  for  knowing  the  characteristics  of  the 
solubility  of  hydrogen  in  palladium;  adsorption  processes  at  the  surface 
of  the  palladium",  which  will  tell  whether  a  certain  impurity  would 
chemisorb  or  physically  adsorb  at  the  surface;  atomistic  diffusion 


t 
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mechanisms  inside  the  solid,  which  deal  with  atomic  dif£u8ipn  through 
different  crystallographic  plane;  and  finally  the  effect  of  dislocation 
in  the  solid  on  the  diffusion  rates.  Some  of  these  discussions  will 
be  qualitative  in  nature,  but  the  author  trusts  that  they  will  serve 
as  a  catalyst  for  future  research. 

B.  Previous  Work 

As  mentioned  previously,  the  transport  rate  of  pure  hydrogen 
in  pure  palladium  produced  numerous  inconclusive  data  mainly  because 
the  investigators  used  static  systems  instead  of  continuous  systems. 
Lombard  and  Eichner  (8)  were  among  the  first  to  use  a  continuous  system 
for  obtaining  data;  however,  their  transport  rates  were  very  or 
relatively  low  owing  to  contamination  of  the  membrane  by  brazing  with 
silver  solder.  Darling  (7)  using  commercial  hydrogen  showed  that,  for 
a  static  system,  the  transport  rate  decreased  rapidly  with  time;  he 
attributed  these  results  to  the  fact  that  impurities  blocked  the 
entrance  of  hydroged  at  the  surface.  When  a  bleed  was  imposed,  the 
transport  rate  did  not  vary  with  time,  thus  lending  evidence  to  his 
thesis . 

Hurlbert  and  Konecny  (5),  using  a  continuous  system,  observed 
that  surface  reaction  is  not  the  rate-limiting  process  for  membranes 
thicker  than  20  microns.  DeRosset  (6)  used  pure  palladium  to  study 
the  transport  rate  of  impure  feed  gas  with  an  upstream  pressure  of 
1-700  psig  and  a  downstream  pressure  of  0-300  psig. 

Commercially,  Engelhard  Industries  (9)  and  Milton  Roy  both 
use  palladium-silver  alloy  as  the  separation  barrier  for  their  small 
scale  hydrogen  purifier.  Union  Carbide  (10)  has  constructed  a  large 


scale  hydrogen  purifier  using  palladium  as  membrane,  but  few  operating 
data  are  given,  A  mixed  feedstock  of  gases  can  be  used  In  all  of 
these  commercial  hydrogen  purifiers. 

^  Pure  palladium  membranes  have  the  disadvantage  of  cracking  and 

distorting  when  used  for  a  long  time,  especially  when  the  membrane  is 
cycled  through  heating  and  cooling  In  an  atmosphere  of  hydrogen. 
However,  Union  Carbide  uses  Its  palladium  membrane  continuously  without 
shuting  dowa,   and  the  life  of  the  palladium  is,  therefore,  prolonged. 
Pure  palladium  usually  results  in  a  higher  transport  rate  at  high 
pressures  and  medium  temperature  than  palladium-silver  alloy,  because 
the  maximum  solubility  of  hydrogen  in  pure  palladium  is  greater  than 

^  palladium-silver  alloy.  Engelhard  Industries  and  Milton  Roy  utilize 

a  palladium- stiver  alloy  in  their  diffuser.   Since  these  units  are 
usually  small  and  are  often  shut  off  and  turn  on,  the  palladium  alloy 
has  an  edge  over  pure  palladium. 

The  bulk  of  the  measurements  of  hydrogen  transport  through 
palladium  are  for  pure  palladium  membranes  and  for  a  downstream 
pressure  of  greater  than  an  atmosphere.  When  divergent  results  were 
obtained  by  the  investigator  for  the  transport  of  hydrogen  through 
palladium  usually  qualitative  explanations  were  given. 

c 

-^  McBride  and  McKlnley  (10)  found  that  when  carbon  monoxide, 

methane,  or  hydrogen  sulfides  were  present  in  the  feedstock,  they 
reacted  with  palladivnn  at  a  temperature  slightly  above  400**C  to  form 

I 

palladium  carbides  and  at  a  lower  temperature  to  form  palladium 
sulfides.  Hurlbert  and  Konecny  (5)  observed  sulfide  formation  at 
■  around  350  C.  Several  investigators  (5,11)  noticed  that  air  will 
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reactivate  a  poisoned  palladium  membrane  when  heated  at  high  temperature; 
this  the  author  believed  to  be  a  removal  of  poison  by  oxidation.  The 
palladium  oxide  formed  in  this  case  can  be  easily  removed  when  the 
membrane  Is  exposed  to  an  atmosphere  o£  hydrogen,  and  this  process 
will  be  examined  thermodynamically  in  a  later  section,  inhere  also 
exists  a  possibility  that  carbon  is  deposited  as  a  result  of  a 
dehydrogenatlon  process,  in  which  the  surface  carbon  may  also  diffuse 
inside  the  metal,  thus  further  poisoning  the  membrane. 

C,  Statement  of  the  Problem 

The  primary  object  of  this  investigation  was  to  study  the 
effect  of  major  impurities  on  the  transport  of  hydrogen  through 
solid  palladium-silver  electrodes.  Hydrogen-oxygen  fuel  cells 
generally  use  porous  material  as  electrodes  and  pure  hydrogen  as  fuel. 
The  use  of  solid  electrodes  has  the  advantage  of  preventing  the  pene* 
tration  of  the  electrolyte  into  the  electrode.  The  use  of  Impure 
hydrogen  as  fuel  would  reduce  the  operating  cost  considerably. 

Experimental  data  were  to  be  gathered  on  the  effect  of  impuri- 
ties such  as  CO^,  CO,  N-,  and  CH,  on  the  transport  rate  with  temperature, 
pressure  and  bleed  rate  as  parameters.  These  impurities  might  be 
expected  to  be  present  in  substantial  quantities  in  the  effluents  from 

ammonia  or  hydrocarbon  reformers.  The  temperature  range  covered 

o       o 
extended  from  300  C  to  460  C  and  pressure  range  from  1069  mm  Hg  to 

2310  mm  Hg. 
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CHAPTER  II 


THEORY 


A.  Thermodynamics  of  the  Pd-H  and  Pd^Ag-H  Systema 

The  rate  of  diffusion  of  hydrogen  through  palladium  is  directly 
proportional  to  the  solubility  of  hydrogen  in  palladium.  Since  the 
equilibrium  isotherms  given information  on  solubility  as  well  as  on 
phase  transitions,  it  is  therefore  important  to  know  the  characteristics 
of  these  equilibrium  isotherms.  Information  on  phase  transitions  Is 
valuable  since  the  diffusion  coefficient  changes  when  a  phase  change 
occurs . 

Sieverts  and  co-workers  (12,13,14)  found  that  hydrogen  adsorbed 

1/2 
in  palladium  according  to  P    law  at  high  temperature  and  low  pressure. 

This  suggested  that  hydrogen  molecules  dissociate  into  atoms  before 

entering  the  crystal  lattice,  indicating  that  hydrogen  possibly  exists 

in  palladium  as  atoms.  However,  Norberg  (15)  and  Isenberg  (16)  proved 

experimentally  and  theoretically  that  hydrogen  exists  in  the  palladium 

lattice  as  quasi- free  protons.  The  electrons  are  attached  more 

closely  to  palladium  atoms  than  to  the  proton.  ; 

Paramagnetic  susceptibility  studies  (17)  show  that  susceptibility 
decreases  with  the  amount  of  hydrogen  absorbed  and  drops  to  zero  when 
there  are  about  0.55  hydrogen  atom  per  palladium  atom.  This  Is 
attributed  to  the  hydrogen  electron  going  into  the  d  shell  of  the 
palladium  atom,  decreasing  the  number  of  positive  holes  and  hence  the 
paramagnetism. 

Gillespie  and  cGalst.aun-  (X8)  obtained  the '.equilibrium 


Isotherms  for  the  Pd-H  system  shown  in  Figure  1.   It  may  be  observed 
•         that  at  a  constant  temperature  and  at  low  pressure,  the  solubility  of 
hydrogen  increases  as  pressure  increases.  Section  A  to  B  in  Figure  1". 

CAt  point  B  the  curve  becomes  horizontal,  indicating  that  the  solubility 
Increases  even  though  the  pressure  remains  constant.  This  increase  in 
solubility  at  constant  pressure  continues  until  a  point  is  reached 
when  an  increase  in  pressure  is  again  necessary  to  increase  the 
solubility,  as  shown  in  Section  CD  in  Figure  1.  The  first  portion  of 
the  isotherm  AB  represents  an  a-phase,  the  horizontal  portion  BC 
indicates  an  a-^  phase  transition,  while  the  last  portion  of  the 
curve,  CD,  represents  the  jB-phase. 

The  length  of  the  a-j3  phase  transition  decreases  as  the 
m-  tanperature  increases  and  finally  disappears  altogether.  The  tempera- 

ture and  pressure  where  it  just  disappears  is  called  the  critical 
point,  and  for  the  palladium-hydrogen  system  it  is  around  20  atmospheres 
and  295°C. 

Nace  and  Aston  (19,20)  studied  the  heat  capacity  of  Pd.H  at 
low  temperatures  ranging  from  16  to  340°K  and  found  that  tetrahedral 
PdH^  is  formed  at  a  temperature  of  less  than  SG^K.  The  PdH,  structure 
is  located  at  the  corner  of  the  palladium  lattice.  At  temperatures 
^  higher  than  50  K  the  PdH  structure  dissociates  into  simpler  hydride, 

and  at  room  temperature  this  model  will  be  Indistinguishable  from  the 
quasi- free  proton  model.  It  can  be  concluded,  therefore,  that  above 
room  temperature,  the  protons  are  distributed  randomly  at  the 
interstitial  sites. 

The  lattice  parameter  of  the  a-phase  is  3.983X;  this  expands  to 
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Figure  1.   Equilibrium  Isotherm 
for  the  Pd-H  Systems 
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4.025A  £or  Che  /3-phase. 

To  Improve  the  mechanical  properties  o£  the  palladium  membrane 
In  an  atmosphere  o£  hydrogen,  palladium  Is  alloyed  with  silver. 
Introduction  of  silver  Into  the  palladium  lattice  changes  the  equili- 
brium Isotherm.  Although  the  equilibrium  Isotherm  for  the  Pd-Ag-H 
systems  Is  similar  In  shape  to  the  Pd-H  system,  the  critical  point 
la  lower  and  the  maximum  solubility  Is  less  than  for  pure  palladium. 
It  was  mentioned  earlier  that  the  electron  of  the  hydrogen  atom  goes 
Into  the  d-band  of  the  palladium  atom.  Since  the  added  silver  atoms 
fill  some  of  the  d-band  vacancies  In  the  palladium  atoms  with  their 
own  electron,  it  will  decrease  the  maximum  solubility  of  hydrogen  In  the 
alloy.  However,  in  the  a-phase  region  at  low  temperature,  hydrogen  is 
more  soluble  in  the  alloy  than  in  pure  palladium.  Figure  2  (39)  shows 
this  characteristic  for  a  25  C  Isotherm.  It  may  be  noted  that  for  this 
particular  Isotherm,  the  equilibrium  pressure  of  hydrogen  for  the 
a-^  phase  transition  diminishes  with  Increase  in  sllyer  content  up  to 
30  per  cent.  Very  few  data  are  available  for  the  Pd-Ag-H  system, 
expecially  at  high  temperature. 
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Figure  2.   Equilibrium  Isotherms 
for  Ag-Pd-H  System  at,  25^0 
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°-   statistical  Thermodynamics  and  Rgulllbrlum  I«n^hpr^,a 

Lacher  (21)  obtained  a  theoretical  equation  for  the  equili- 
brium Isotherms  by  statistical  thermodynamics  based  on  an  Interstitial 
hydrogen  model.  Llbowltz  (22)  arrived  at  a  similar  equation  using  a 
9  vacancy  model.  However,  the  Lacher  model  Is  a  better  model  since  It 

Is  based  on  a  realistic  quasi-free  proton  model  which  has  been 
confirmed  by  several  investigators. 

Here  an  equation  will  be  derived  which  is  a  combination  of  the 
models  used  by  Lacher  and  Fowler  (21,37).  The  assumptions  made  in 
the  derivation  are: 

(1)  There  is  a  fixed  number  of  interstitial  sites  N 

'  s* 

where  the  hydrogen  atom  can  be  absorbed. 
^  (2)  The  energy  of  interaction  between  a  pair  of  nearest 

neighbors  of  hydrogen  is  equal  for  all  pairs. 

Let  the  number  of  pairs  of  neighboring  hydrogen  atoms  in  a 
particular  configuration  be  Vl^,   and  let  2l(^^e^^  be  its  interaction 


z 


energy.  Then  from  statistical  thermodynamics  (37),  the  expression 
for  the  partition  function  for  this  system  will  be: 

2e 

-N    -sa 

r-1  HH  zkT  r      -,N„ 

^"l^CWexp        [a^(T)j«  (1) 

where   Q  -  partition  function 

2^*'h'*'hH^  "  ^°^^^   """**"  of  distinguishable  configurations  of  N  atoms 

H 

z  ">  number  of  nearest  neighbors 
Njj  »  number  of  hydrogen  atoms  in  the  lattice. 
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Euj,  ■  total  Interaction  energy  at  saturation  per  atom 
k  B  Boltzmann  Constant 
T  B  temperature 

•  s^Yi^T!)  "  the  partition  function  for  the  Internal  degree  of 

freedom  of  an  absorbed  hydrogen  atom  referred  to  the 
ground  state  at  Infinite  separation 

la 

Defining  H       by  Equation  (2) 

exp  2^  8(N„.Njjg) 

r-1  -N  2c  /zkT 

-  2^g(Ng,Njjjj)exp  ™  ™  (2) 

^  Substituting  Equation  (2)  into  (1)  the  partition  function  will  be 

in  the  form 


Q  - 


Nhh 


*''»'  ^  N  1 


>  ^%'V   "   N'(/-N)!  (*) 


N  tl    S     H 


therefore.  Equation  (3)  becomes 


Q  -  r  a_(T)  1  "exp*  ™  ^»»  '  ^s'  (5) 

L        J  V<«s-V= 

From  Statistical  thermodynamics  the  Heloholtz  free  energy, 
F,  can  be  V;.' It  ten  ast 

F  »  -  kTlnQ  (6) 

Substituting  the  value  of  Q  from  Equation  (5)  to  equation 


• 
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12 
(6)  and  using  Sterling's  approximation 


F/kT  =  -NglnNg  +  NjjInNjj  +  (H^  -  Njj)ln(N^  -  Njj)   -  Nh^°«H^^^ 


^  VW^^^  <7) 


The  equilibrium  or  average  value  of  N,_,  can  be  expressed 

nn 


from  statistical  thermodynamics  as: 


NffH 


N.„.  =  fllSi 

"^'HH'""'HH' 


HH     TT- — :F-2r-7^kT—  <8) 


2_^g(Njj.Nj^j)exp 

Each  occupied  site  has  z  neighbors,  and  with  perfectly  random 

arrangements  each  site  has  a  probability  rr-  of  being  occupied. 

"s  N 

Therefore  the  average  number  of  neighbors  of  any  given  N„  Is  z  rp^.  The 

H      N 
S 

total  equilibrium  number  of  pairs  of  neighbors  will  be: 

"kh=(1''h)(==f)447   ■        O 

But  N   =  N   for  random  distributions  (38);  therefore 
Equation  (9)  becomes 

It  =  "^s^^^s  -^  V^^H  ^  <«s  -  V^"<»'s  -  V  ■  V"^h(^> 

-Fkr^  (10) 

s 
The  chemical  potential  can  be  written  as: 
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296, 


No 
where  ^  ■"  ^  the  fraction  of  sites  occupied  but 


•e„/kT 
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^  aflCT)  -  exp  "   qjj(T)p  (12) 

where 

(-Cjj)  is  the  energy  required  to  remove  an  absorbed  H  atom,  far 
removed  from  any  other  absorbed  H-atoms,  from  its  lowest  state  in 
the  metal  to  rest  at  infinite  dispersion. 

qjj(T)  is  the  partition  function  for  vibration  of  absorbed 
H  atom  relative  to  the  lattice. 

p  is  the  spin  weight  of  the  proton. 
^  Defining  \„  as  the  absolute  activity  and  using  Equation  (12) 


V°>  ■  -^"^  -  rh  Trx).  (») 


The  absolute  activity  in  the  gas  phase  consisting  entirely  of 
diatomic  molecules  of  hydrogen  at  a  pressure  f>   is: 

,-,  .1/2   -1/2V"  r(2„™kl)3'2  M^Wl"" 

(14) 

where  e^  ■>  energy  required  to  dissociate  a  hydrogen  molecule  from 
its  lowest  state  into  two  free  H-atoms  each  in  its 
lowest  state. 

IL,  'o  atomic  mass  of  hydrogen 

h  «  Planck  Constant 

M-  a  isoment  of  iniertla  of  the  hydrogen  molecule. 
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Equating  Equations  (13)  and   (14) 


/     0      \2  2(€„  +  2%^  +  ^  Jk/T       (27l2M„kT)3/2 


87r^MjkT 
'^       2  r     ,„a  2  <15) 


2h'  fq„(T)3 
Let  P^  2   )  '>•  defined  by   Equation  (16) 


(^y     {q„(T)3^     ~V'  ir~     '''' 


3/2  „  2„  ,,  ^^^H  ^  ^HH  -^  K> 


kT (27l2M^kT)"^  "  87r"M.|.kT  2 ^ 

1 


2(e„  +  e      +  ^„,)/kT     ■• 
-  K(T)exp       "         HH       ^  D  (ig) 

where 

,_  (27i2M„kT)^^^  STT^M^kT 

^q  (T)J  ^  h^  2h2 


Therefore  Equation  (15)  will  be  simplified  to: 


2       (20-   l)2e_/kT 
^"Mt   ){rf-5   )exp  ««  (18) 


If  ^Tni^O"  there  will  be  a  critical  temperature  given  by: 

T  "  critical  temperature 
It  is  not  possible  to  determine  theoretically  all  of  the 
parameters  in  Equation  (15),  therefore  experimental  data  must  be 
obtained  to  evaluate  these  constants.  To  do  this  it  is  convenient 
to  convert  Equations  (16),  (18)  and  (19)  to  a  molal  basis. 
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(I  j-K'(T)ejcp   H   Tffl   2  D  ^^o) 


0   v2   (2©-l)2E  /RT 
P-Pf  f  )f  t-~  Uxp      ^"  (21) 


"V'^'^c  "2  ,   (22) 

Th«  constants  K'  (T)  and  E^  +  E^  +  ^  E^  are  obtained  from  the 
plot  of  InP  [^  \  versus  -.  While  Ey^^  Is  obtained  from  the  critical 
temperature  by  using  Equation  (22). 

Lacher  (21)  obtained  the  following  equation  similar  to  equation 
(20) 

l°8^n»i(i)'=  7.4826-^221^  (23) 

He  also  found  that  T^  -  568°K.  so  that  Ejjjj--2(R)(568),  Taking  the 
log  of  Equation  (21)  and  substituting  the  corresponding  value  of 
E^  one  obtains 

On  inserting  Lacher 's  value  for  log  Pf  ^  j  Equation  (24)  becomes 

logP^  -  7.4826  -  -iSIIJl^  +  21og  -A-  .  986.7(26  -  1) 
mm  X         1  "  D       T 

-  7.A826  f  21,g -^  ■'"•';  ""-AS      (25) 

Equation  (25)  is  plotted  in  Figure  1  along  with  typical 
experimental  data.  It  can  be  seen  that  Equation  (25)  correlates  very 
well  with  the  data. 

Equation  (25)  is  obtained  by  using 'pure  palladium  data, 
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therefore  it  cannot  be  used  for  the  palladium- a liver  alloy.  However, 
the  few  experimental  data  available  for  the  palladium-silver  alloy 
suggest  a  similar  pattern.  Since  there  is  a  lack  of  data  for  the 
whole  range  of  the  equilibrium  isotherms  for  the  palladium-silver- 
hydrogen  system,  an  extrapolation  method  will  be  used  to  predict   « 
these  equilibrium  isotherms. 

The  equilibrium  isotherms  for  the  Pd-Ag-H  system  at  25°C 
for  several  silver  contents  were  shown  in  Figure  2.   It  can  be  seen 
that  the  equilibrium  pressure  for  the  a-A  phase  transition  decreases 
in  proportion  to  the  silver  content.  This  phenomenon  can  be  used 
to  predict  some  of  the  constants  in  Equation  (21)  for  the  Pd-Ag-H 
system. 

In  Figure  3  a  plot  of  logP^j^^j^  (jj  versus  the  silver  content 
of  the  alloy  obtained  from  Figure  2  at  25°C  yields  the  following 
relationship 

1°8%,  "  1-2165  -  6.52X^g    ,.  (26) 

where  X   is  weight  fraction  of  silver. 

Assuming  that  this  relation  applies  also  to  other  temperatures, 


so  that 

a 


Ij  %  6.52  when  T  =  298°C  (27) 

a  o  1942 


and 

l02 

"mm  \ZJ        '.-'--         ,j.       «.^g 


log  ~P{~]-  1.2165  -  4^  X. 


(0- 

iMerting  this  value  for  log  P^^  fj\  in  Equation  (23)  we  obtain 
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log  P^  (I)  -  7.4826  -  M2^  -  ^  X. 


'tad  V2;    '••*°'°    ^     T  \g 


^  7.4826 5 ^S-  (29) 

^  For  the  prediction  of  t^   for  the  Pd-Ag-H  system,  a  plot  of 

T^  versus  X   Is  made.  Only  two  points  are  known  for  this  plot;  one 
of  the  points  Is  obtained  from  the  data  on  the  30%  Agof  Makrides  (23), 
and  the  other  Is  for  pure  palladium.  Assuming  a  straight  line 
relationship  Figure  4  was  constructed. 

One  conclusion  which  can  be  shown  from  Figure  4  Is  that  \, 
decreases  as  the  amount  of  silver  Increases.  This  can  be  explained  ' 
by  the  fact  that  the  total  Interaction  energy  at  saturation  per  mole 

^  will  decrease  because  the  solubility  of  hydrogen  Is  decreased  by  the 

addition  of  silver.  Makrides  (23)  used  the  ratio  of  maximum  solubility 

of  hydrogen  In  alloy  to  that  In  palladium  to  correct  the  value  of 

E   In  the  alloy. 
HH 

The  value  of  T  for  the  25%  silver  alloy  Is  about  350°K: 
c  ' 

therefore  the  equilibrium  Isotherm  expression  for  the  25%  silver  alloy 
will  be: 


^°8^tem  -  7-^826  +  21og  ^  -  1635.52  +  12180     ^^^^ 

C.  Transport  of  Hydrogen  Through  Palladium 
The  transport  of  hydrogen  through  palladium  can  be  attributed 

to  the  following  series  of  processes: 

(1)  The  hydrogen  gas  molecule  is  adsorbed  and  dissociates 

Into  hydrogen  atoms  at  the  upstream  surface. 
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Figure  4,   A  Plot  of  Critical  Temperature 
Versus  Silver  Content 
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(2)  The  adsorbed  hydrogen  atom  enters  Into  the  crystal 
lattice  and  becomes  a  proton  and  a  quasi-free  electron. 

(3)  The  proton  diffuses  through  the  palladium  lattice. 

(4)  At  the  downstream  surface  the  proton  recombines  with 
^             an  electron  to  form  a  hydrogen  atom. 

(5)  TWO  hydrogen  atoms  combine  to  form  a  hydrogen  molecule. 

(6)  Hydrogen  molecules  desorb  from  the  downstream  surface. 
Each  of  the  steps  above  may  be  the  rate- control  ling  one;  it 

is  therefore  important  to  discuss  them  separately. 

1.  The^  Adsorption  and  Dissociation  of  ........   .....  „^^^^^ 

The  hydrogen  molecule  before  diffusing  through  the  metal  must 
Q  first  adsorb  in  the  metal  surface  and  subsequently  dissociate  into 

two  hydrogen  atoms.  Adsorption  processes  can  be  broadly  divided 
into  two  types  depending  on  the  amount  of  heat  evolved  during 
adsorption;  these  are  (1)  physical  adsorption,  (2)  chemisorption. 

The  heat  of  chemisorption  is  generally  higher  than  the  heat 
Of  Physical  adsorption,  since  in  chemisorption  chemical  bonds  appear 
to  be  formed  between  the  gas  and  metal  atom.   In  physical  adsorption 
van  de  Waals  forces  are  the  forces  acting  between  molecules.and 
^  these  are  generally  small;  thus  the  heat  evolved  is  generally  small. 

Electrochemical  determination  of  the  heat  of  adsorption  of  hydrogen 
or  finely  divided  palladium  shows  a  heat  of  adsorption  of  27.5  kcal/mole 
(24),  which  is  in  the  range  of  chemisorption. 

Chemisorption  is  usually  the  first  step  toward  the  dissocia- 
tion  Of  the  hydrogen  molecule  on  metal  surfaces;  it  is  found  to  be 
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true  especially  In  heterogeneous  catalytic  chemical  reactions. 

The  conditions  in  which  surface  effects  determine  the  rate 
of  diffusion  are  the  following; 

(1)  low  pressure  and  low  temperature, 
'I*.  (2)  very  thin  membrane, 

(3)  blocking  and  poisoning  of  the  surface  by  impurities, 

(A)  physical  condition  of  the  surface. 

Hydrogen  molecules  are  adsorbed  and  dissociated  on  certain 
active  sites  in  the  metal  surface,  the  number  of  molecules  adsorbed 
and  dissociated  being  proportional  to  the  product  of  the  number  of 
molecules  striking  the  surface  and  the  probability  of  adsorption 
and  dissociation.  This  condition  therefore  requires  relatively  high 
Q  pressure  and  temperature.  At  low  pressure  the  number  of  molecules 

striking  is  greatly  reduced,  resulting  in  a  reduction  of  the  trans- 
port rate.  The  rate  of  dissociation  increases  exponentially  with 
temperature,  and  proportionally  with  square  root  of  pressure.  This 
means  that  temperature  should  produce  more  profound  effect  than 
pressure. 

The  surface  effect  will  also  be  the  rate-controlling  factor 
if  the  thickness  of  the  metal  membrane  is  greatly  reduced  so  that 
^  the  resistance  to  flow  offered  by  the  metal  membrane  is  negligible 

compared  to  that  of  the  surface.  This  means  that  the  metal  membrane 
will  diffuse  all  the  gases  that  the  surface  has  dissociated.   Hurlbert 
and  Kbnecny  (5)  found  that  at  upstream  pressures  of  1  to  7  atmospheres 
and  temperature  around  350°C,  surface  effects  are  rate-limiting  for 
membranes  thinner  than  20  micron.  This  limiting  thickness,  of  course, 
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will  vary  for  different  temperatures and  pressures. 

The  blocking  and  poisoning  of  the  surface  by  Impurities  Is 
caused  by  occupation  of  the  active  sites  by  foreign  atoms  (25,26). 
The  following  kinds  of  substances  can  cause  poisoning: 

(1)  impurities  in  the  feed  gas, 

(2)  vaporized  metal, 

(3)  impurities  in  the  metal, 

(4)  metal  oxide  formation. 

Impurities  in  the  feed  gas,  such  as  hydrocarbon  vapor  and 
hydrogen  sulfide,  decrease  the  transport  rate  greatly  at  temperatures 
In  the  range  400°C  to  600°C  (25).  This  was  explained  by  McBrlde  (10) 
as  caused  by  the  chemical  reaction  of  these  two  gases  with  palladium 
forming  palladium  sulfide  or  carbide.  The  formation  of  palladivmi 
sulfide  is  instantaneous  and  irreversible.  The  formation  of  palladium 
carbide  may  involve  a  more  complicated  reaction,  and  may  be  a  very 
slow  process;  furthermore  the  carbon  atoms  at  the  surface  may  diffuse 
into  the  palladium  lattice,  thereby  further  decreasing  the  transport 
rate. 

Oxygen  and  water  react  with  palladium  to  form  palladium  oxide. 
If  an  atmosphere  of  hydrogen  gas  is  introduced,  the  following  reaction 
will  occur 

PdO  +  H^  ^  Pd  +  H2O  (31) 

From  thermodynamics 

log  K  =  -AF°^/2.3RT  (32) 

where   K  »  the  equilibrium  constant  for  a  reaction. 

AF°*  -  the  free  energy  of  formation  at  298°k  in  cal/mole. 
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For  the  reaction  above 

AF°  =>  AF  (product)  -  AF  (reactant) 

-  -54.64  -  (-13.45)  =  -41.19  kcal. 
Substituting  the  value  of  AF  and  a  temperature  of  623*'k 
to  equation  (32) 

logK  =14.5 

K  =  3.2  X  10^^ 

\^  14 

~-  =»  3.2  X  10^^  (33) 


■«2 


Therefore 


%  •  '  ■=  1°''  \o 


Therefore  if  the  vapor  pressure  of  water  is  around  23  mm  Hg, 

the  corresponding  equilibrium  pressure  for  hydrogen  will  be 

-13 
1.2  X  10   mm  Hg;  since  the  pressure  of  hydrogen  is  usually  above 

760  mm  Hg,  it  can  be  concluded  that  all  of  the  PdO  will  be 

converted  to  palladium. 

Similarly  for  the  silver,  the  reaction  will  be 


Ag20  +  Hg  ^  2Ag  +  HgO  (34) 

From  this  it  can  be  concluded  that  no  Ag.O  will  exist  when  a 
hydrogen  pressure  of  more  than  760  mm  Is  present. 

Many  Investigators  have  found  that  an  inactive  palladium 
membrane. can  be  reactivated  by  heating  in  an  atmosphere  of  air; 
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however,  this  method  does  not  produce  a  lasting  effect.  Recently, 
the  use  of  low  energy  electron  diffraction  (27)  has  led  to  some 
understanding  of  the  situation  when  certain  gases  are  adsorbed  on 
the  metal  surface.  It  was  found  that  the  surface  atoms  rearrange 
when  a  gas  is  adsorbed.  The  adsorption  of  oxygen  on  nickel,  for 
example,  causes  rearrangement  of  all  three  crystallographic  planes 
of  nickel  (100,  110,  111  surfaces).  Since  palladium  Is  similar  to 
nickel  in  crystalline  structure,  i.e.,  both  are  face  centered  cubic 
and  both  are  transition  metals,  it  seems  reasonable  to  assume  that 
adsorption  of  oxygen  causes  surface  rearrangement.   In  these 
circumstances  the  transport  rate  will  be  relatively  high  at  the 
start;  then  when  the  surface  rearranges  to  the  equilibrium  configura- 
tion, the  transport  rate  will  decrease.  This  was  observed  by 
Makrides  (23),  but  he  offered  a  different  explanation,  namely  that 
the  decrease  in  transport  rate  was  caused  by  impurities. 

When  the  diffusion  cell  is  operated  at  high  temperatures  some 
of  the  metal  may  be  vaporized  and  is  deposited  on  the  surface  of  the 
membrane.   This  would  be  expected  to  decrease  the  transport  rate. 
Such  metals  as  lead  and  iron  in  the  form  of  FeCl_  can  poison  the 
palladium  membrane.  Surface  diffusion  of  foreign  atoms  from  the 
vicinity  of  the  palladium  membrane,  such  as  copper  gasket,  also  slows 
down  the  transport  rate. 

Assuming  that  there  is  no  effect  exerted  by  neighboring 
molecules,  the  time  in  which  a  molecule  stays  at  the  surface  can  be 
predicted  from  the  heats  of  adsorption  and  the  physical  properties 
of  the  metal.  This  residence  time  of  the  molecule  at  the  surface  can 
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be  used  as  a  measure  of  the  poisoning  effect  of  certain  molecules. 

According  to  De  Boer  (28)  the  adsorption  time  of  a  molecule 
In  the  surface  Is  given  by: 

All 

•^  -  ^o^'^P  RT  (35) 

2/3 
where   T^  »  4.75  x  lO"^^  ^  (36) 

B 

T  =  time  of  adsorption 

M  =  mean  molecular  weight  per  constituent  atom 

3 
V  =>  molar  volume  In  cm 

T  =  melting  temperature  K 

_  o 

T  =»  temperature  K 

R  =  gas  constant,  cal/mole**K 
AH  =>  heats  of  adsorption,  cal 

The  residence  time  for  hydrogen  on  palladium  at  room  tempera- 

-2 
ture  Is  about  10  sees.  If  a  gas  has  a  much  larger  adsorption  time 

than  hydrogen,  then  the  gas  can  be  said  to  be  a  poison. 

Trapnell  (29)  Indicated  that  unsaturated  hydrocarbons,  such  as 

ethylene,  chemlsorbed  on  the  transition  metals  according  to  the  equation: 

CH  "  CH 
I    / 
4Pd  +  CgH^-*  Pd  -  Pd  +  2PdH  (37) 

The  heat  of  adsorption  In  this  case  can  be  estimated  by  using 
Pauling's  (30)  equation  for  single  bond  energies. 

E(X  -  Y)  »  I  E(X  -  X)  +  E  (Y  -  Y)  +  2.306(/ijj  -  Py)^       (38) 
where  E(X  -  Y)  is  the  bond  energy  between  X  and  Y 
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^  Is  the  electronegativity. 

The  initial  heat  of  chemisorption  q  can  now  be  written  as: 

q^  "   2E(Pd-Pd)  +  E(H-H)  +  E(C-C)  -  2E(C-H) 

Substituting  the  values 

E(Pd-Pd)   -  15.5 
E(H-H)        -  104.2 
E(C-C)        «  83.1 
E(C-H)        =  98.8 

/4d-C       =°-°^ 

^  4d-H         ='°-01 

one  obtains  a  value  q^  of  22.5  kcal  and  an  adsorption  time  of 

3 
4  X  10  sees.  Similarly,  this  approach  yields  an  adsorption  time  for 

hydrogen  of  4  x  lO'  sees.  Therefore  it  can  be  concluded, that  for 

an  equal  mole  fraction  of  hydrogen  and  ethylene  at  room  temperature,  the 

ethylene  will  be  adsorbed  more  than  hydrogen.  However,  the  acetylinic 

complex  formed  in  Equation  (37)  is  slowly  removed  by  gaseous  hydrogen, 

presumably  as  ethylene: 

A  CH  "  CH 

Pd  -  Pd  +  H2  -^  2Pd  +  C2H^  (40) 

or 

H    H 

I    I 

Pd  -  Pd  +  C  H  -H.  2Pd  +  C  H  (41) 
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2.  The  Diffusion  of  Hydrogen  In  the  Palladium  Lattice 

There  are  several  theories  which  attempt  to  explain  the  absorption 
of  hydrogen  In  palladium  and  the  mechanism  of  Its  diffusion  through  the 
palladium  lattice.  Some  of  these  theories  are: 

(1)  A  chemical  compound  Is  formed  when  hydrogen  Is  absorbed 
in  the  palladium  lattice. 

(2)  Hydrogen  Is  rift  occluded  In  the  palladium  lattice. 

(3)  Hydrogen  forms  an  Interstitial  solution  either  as  atoms 
or  protons  In  the  palladium  lattice. 

Simons  and  Ham  (31)  considered  the  diffusion  of  gases  through 
metal  as  pure  chemical  phenomenon,  I.e.,  a  chemical  compound  is  formed 
between  the  gas  and  the  metal.   This  theory  would  suggest  that  diffusion 
of  hydrogen  Is  accomplished  by  the  kinetic  dissociation  of  the  compound, 
and,  since  a  concentration  gradient  exists  in  the  metal  membrane,  the 
movement  of  hydrogen  will  be  toward  the  low  pressure  side.  This  theory 
may  be  valid  if- the  bond  between  the  metal  and  gas  is  sufficiently 

^  strong  to  consider  it  a  compound.   However,  since  hydrogen  diffuses 

through  palladium  quite  rapidly  at  elevated  temperatures,  it  appears 
that  this  mechanism  is  of  minor  importance.  Nace  (19,20)  has  suggested 
that  chemical  bonds  are  formed  between  the  hydrogen  and  palladium.  The 

(  )  hydrogen  diffuses  through  palladium  in  this  case  is  assumed  to  occur  as 

the  movement  of  hydrogen  atoms  by  the  following  mechanism:   (1)  rotation 
about  a  parent  palladium  atom;  (2)  bonding  to,  and  subsequent  rotation 
about  a  neighboring  atom.  When  the  temperature  is  raised  to  room 
temperature,  the  bonds  are  broken,  and  the  hydrogen  atoms  are  randomly 
distributed  so  that  an  interstitial  solution  is  formed. 
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Smith  (32)  believes  that  hydrogen  diffuses  through  palladium 
via  rifts,  a  rtft  being  a  deformed  or  dislocated  portion  of  the  metal 
crystal  possessing  a  great  Ionization  power  because  of  Imbalanced 
forces  In  the  stressed  metal.  This  theory  was  concluded  to  explain 
the  greater  solubilities  and  diffusion  rates  observed  when  the  metal  is 
deformed  by  cold  working.   It  is  believed  that  the  formation  of  rifts, 
as  shown  by  line-broadening  in  the  X-ray  diffraction  spectra,  causes 
hydrogen  to  be  occluded  in  the  metal.  Kydonlus  (33)  found  that  sonic 
vibration  Increased  the  rate  of  mass  transfer  through  palladium  by 
producing  dislocations.   It  would  be  hazardous  to  assume  that  the 
mechanism  of  diffusion  is  mainly  through  rifts  in  dislocation.  Because 
If  this  Is  true,  other  metals  such  as  nickel  and  copper  can  be  made  to 
have  the  same  transport  rate  as  that  of  palladium;  this  has  never  been 
accomplished  experimentally.  The  probable  reason  for  the  Increase  in 
the  absorption  or  diffusion  rate  in  stressed  metal  is  that  In  stressed 
metal  the  dislocations  serve  as  a  reservoir  for  the  hydrogen  and  at  the 
surface  dislocation  Increases  the  surface  activity.  However,  the  majority 
of  hydrogen  appears  to  diffuse  inter s tit la lly. 

The  diffusion  of  hydrogen  in  the  palladium  lattice  can  be 
described  from  an  atomistic  point  of  view.   It  was  shown  earlier  that 
the  proton  occupies  certain  interstitial  sites  in  the  palladium  lattice, 
the  most  probable  sites  being  those  shown  in  Figure  5.  These  sites 
are  located  at  the  edges  and  at  the  center  of  the  face-centered  cubic 
lattice. 

Diffusion  appears  to  involve  hydrogen  jumping  from  one  inter- 
stitial site  to  another.  From  these  considerations  an  equation  can  be 
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Figure  5.   Interstitial  Sites  for  a  Face  Center  Cubic  Lattice. 
0  Indicates  Metal  Atom,  X  Indicates  the  Interstitial  Position 
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derived  to  describe  these  processes  (40,41).   Consider  the  interstitial 

position  of  a  face-centered  cubic  lattice  as  shown  in  Figure  6.   Three 

atomic  planes  are  drawn  in  this  figure  separated  by  a  distance  of  X. 

Assume  that  there  exists  a  concentration  gradient  of  hydrogen  along  the 

X  axis  which  is  perpendicular  to  the  (001)  atomic  plane.   A  hydrogen 

in  an  interstitial  position  may  jump  in  forward,  backward,  or  in  a 

direction  perpendicular  to  the  x-axis.   Let  P  be  the  probability  for 

a  given  interstitial  hydrogen  to  make  a  jump  per  second  and  f  P  be  the 

r  r 

probability  for  a  jump  per  second  in  the  forward  direction.   Also 

assume  the  probability  for  a  forward  and  backward  jump  to  be  equal.  Let 

2 

the  number  of  diffusing  particles  per  cm  on  the  plane  located  at  x  at 
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Figure  6.     Atomic  Plane  of  Diffusion 
in  the  Metal  Lattice 
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the  Instant  t  be  represented  by  n(x).  Using  Taylor  series  expansion 
one  therefore  gets 


n(x  +  X)  "  n(x)  + 


(^>^i(&y^ 


n(x-X)-nU).(|sy4(^^y 


(42) 


(43) 


lAere   n(x  +  \)  •=■  the  number  of  hydrogen  per  cm  at  plane  X  +  X 

n(x  -  X)  ■  the  number  of  hydrogen  per  cm  at  plane  x  -  X 
at  time  t 

The  increase  of  the  number  of  particles  located  at  plane  x  at  instant 
t  +  6t  will  be: 
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6n(ic)  =.  r  f^P^n(x  +  X)  +  fyP^nCx  -  X)  -  2f^P^n(5c)l  6t  (44) 

2 
6n(x)  -  f,Pj6t  ^  X^  ■  (45) 

from  whence  we  have 

|a  -  f  P  X^  ^  (46) 

Comparing  Equation  (46)  with  Flck's  second  law  we  see  that 

D  -  f^P^X^  (47) 

where   D  =  diffusion  coefficient 

f  =»  a  geometric  factor  and  for  a  f.ti.C.  lattice  is  equal 
to  1/3  when  the  atomic  plane  is  001  plane. 

The  probability  P  is  given  by  (36) 

P  =¥^  (A8) 

where   Q+  =>  partition  function  of  the  activated  state 
Q+  3  partition  function  of  the  normal  state. 

Further  simplification  of  equation  (48)  yields 

^  "  T  (^^  -  ^'^  "kT  j  exp  (49) 

where      v  =>  frequency  of  vibration  of  a  linear  oscillator 

AF+  ■=»  free  energy  of  activation. 
For  the  case  where  hl«kT   Equation  (49)  reduces  to 

_         „       -AF+/RT     ,.         -AF/RT  ,,„. 

P^  =  Vexp  «"  V  exp  (50) 

-AF/RT 
D  -  f^vx  exp  (51) 
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For  a  face- centered  cubic  lattice  with  lattice  constant  a, 
and  using  relationship  AF  »  AE  -  IAS  Equation  (51)  becomes 
\2 


3''ft) 


^     AS/R         -AE/RT    -  1         2       AS/R       -M/RT 
exp  exp  ■  12  ^^^  exp"^'   exp  "^''^^ 


».«P-^'"  (52) 


where   AS  =»  entropy  of  activation 

AE  =»  energy  of  activation 

,,12   AS/R 

D^  -  ^  va  exp^"*  (53) 

Calculating  D    with  estimated  values  of  AS,  v  and  a 

V  =■  1.75  X  10^"'  sec."^ 

-8 
a  =  3.9  X  10       cm 


G  AS  -  P(AE/y  -  0.18  C^)   -  0.46 


leads  to 

\-W-  10^3  X  (3.9  X  10-«)2  exp  '-^'^^''^  -  1.75  x  lO-\i/sec. 

(54) 
If  the  hydrogen  is  diffusing  normally  to  001,  010  or  100  plane, 

the  diffusion  constant  would  be  the  same;  however,  if  the  hydrogen  is 
diffusing  normally  to  110  plane,  the  values  of  f  and  X  would  change, 
indicating  that  solids  are,  in  general,  anisotropic.   Initial  calcula- 
tions for  diffusion  through  the  110  plane  Indicate  a  higher  diffusion 
coefficient.  From  these  equations  It  can  also  be  predicted  that  an 
expanded  lattice,  such  as  a  p-phase,  would  have  higher  diffusion 
coefficient  provided  that  no  covalent  bond  la  formed  between  the  hydrogen 
and  the  palladium  atom.  The  diffusion  coefficient  in  Ag-Pd  alloys  (34) 
was  found  to  decrease  with  increasing  silver  content.  This  is 


c 
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anticipated  since  silver  Is  a  poor  hydrogen  dl££user.  One  problem 
that  may  arise  In  the  Pd-Ag  alloy  Is  the  order-disorder  phenomenon 
of  Pd  and  Ag  atom  in  the  lattice.  As  discussed  by  Krivoglaz  and  Smlrnov 

(35>  ,  the  order-disorder  of  the  alloy  will  affect  the  value  of 
energy  of  activation.  This  is  further  complicated  by  the  fact  that 
metal  crystals  are  not  perfect;  therefore,  the  diffusion  would  be  in 
different  crystallographic  planes,  which  were  shown  above  to  possess 
different  diffusion  coefficients. 

3.  Macroscopic  Diffusion  Equation  in  Palladium 

The  macroscopic  equation  for  diffusion  of  hydrogen  in  palladium 
can  be  written  as: 

0  sr  ■"  ^'\  =  \  <55) 

where   C,  =  concentration  of  A 

A 

N,  »  molar  flux  of  A 

A 

R.  =  chemical  reaction  rate  of  A 
A 

V  "*  del  operator 

This  can  be  written  in  rectangular  coordinate  as 

he.        aw..    ^N,    9N, 

5r--^^^^-#  =  -A  (36) 

"^"^^   ^Ax'^-IT-'V  (57) 

V  =•  molar  average  velocity 

»  a.v./2C. 

n  1  i  n  1 

In  the  diffusion  of  hydrogen  through  solid,  v*  will  be  near  zero,  and 
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since  there  is  no  chemical  reaction,  R.  »  0.  In  this  Investigation, 
the  flow  la  In  a  single  direction  only,  say  the  x-dlrectlon;  therefore 
Equation  (56)  becomes 

(  '^'■■^'''  (58) 

>>"'  N.   -  '"^A  (59) 


v^- 


O 


therefore 


'hi 

^*'ax 

Ax 

hx 

dt 

'h(^) 

le  so] 

Lid  to  be  Iso 

3C^        Td\ 

(60) 


(61) 

Let  US  assume  that  the  upstream  surface  reaction  Is  very  fast 
so  that  equilibrium  condition  prevails  between  the  surface  and  hydrogen 
gas.  The  corresponding  boundary  conditions  and  differential  equation 
will  be: 

(62) 
(63) 
(64) 

(65) 

(66) 


(67) 
(68) 
(69) 


B.C.I 

C  =  G^. 

X  =  0, 

t>0 

B.C. 2 

C  =  0, 

X  =  L, 

t>0 

B.C. 3 

C  =  0, 

0<x^, 

t  =  0 

0<1<I'. 

t>0 

Let  V  -  C/C^ 

then 

B.C. 4 

v  -  1, 

X  «=  0, 

t>0 

B.C. 5 

V  =  0, 

x  =  L, 

t>0 

B.C. 6 

V  =>  0, 

Q«L, 

t   =  0 
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^  -  -^  00<L,  t>0  (70) 

ox 

Taking  the  Laplace  transform  of  Eqijatlon  (70)  with  respect  to  t 

2 
(     ,  8V(x,8)  -  v(x.O)  -  ^  ^<^'°>  -  0  (71) 

'■"•  dx 

where  J^ 

V  a     /  vexp"^*^dt  (72) 

o 

since  v(x,0)  •»  0  (73) 

Integrating  Equation  (71) 

V  -  C^exp"*^*  +  C2exp"'l*  (74) 

where  a  =  I   r   I  (75) 


(t 


I'-]' 


Applying  B.C.   (4)  and  (5)  to  eliminate  constants  C.  and  C, 
8(exp^     -   exp       )    L  J 


„  8lnh(L-x)q 

a  slnhqL  ^'°' 

7+='-*  St 

V  =  JL      r  exp     slnh(L-x)q   ^ 

^       27ri     /  s  SlnhqL  "     ._  ^''> 


lOO 


7-lc 
^  The  Inverse  Laplace  transform  of  Equation  (77)  with  poles  of  s  »  0, 

I? 


and  s  o  "2       .  n  "  1,  2,   3  Is 


,^2  2^ 
00  -Pn  Y 

"•*"  ""  -^  '  '  Cos  nrr 


V  -  ilS  +  2    y   exp     /-         sinim^ 


°"1  «  « 

2  2 

^  _         -Dn  TT  t 

L-x 

L 


xt~s»l 


t 
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The  diffusion  flux  across  plane  L  at  any  time  t  will  be: 
-DdCj 


N 

00        -Dn  ~ 


'x-L    dx|^.j 


"^i       2DC,  f    „  =^^ 
-^+-^  2^(-l)"exp  ^  (79) 

n-l 
The  diffusion  flux  at  Infinite  time  will  be 

DC 

N      =  ■—  (80) 

t  — ♦•  00    L 

This  equation  represents  the  steady-state  value  of  the 

diffusion  flux.  Therefore  the  fraction  of  the  steady-state  value 

reached  at  time  t,  SST,  will  be  the  ratio  between  N  j^  and  N^  _^  ^  . 

Thus  the  value  of  SST  becomes 

__        -Dn  TT  t 

SST  "  ;^^=^  -  1  +  2  )  (-l)"exp   ^  (81) 

n^l 

Letting  T  =•  -^^^j  the  value  of  SST  for  different  values  of  T 

Is  plotted  In  Figure  7. 

Figure  8  can  be  used  to  find  the  time  to  reach  907.  of  the 

-6  2 
steady-state  value.  For  example,  if  L  =  0.01  cm,  D  =  10  cm  /sec, 

the  time  t  to  reach  90%  of  equilibrium  value  will  be  (with  T  »  3  when 

f~)  SST  =  0.90  from  Figure  7): 

^„3^(10-^)(3.14)2t  •   ^32) 

(0.01) 

,  and  t  a  30  sees. 

In  the  case  where  there  is  a  surface  effect,  the  equation 

will  be  of  the  form 
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^•^•^  5i"J^^°  -  ^i>    "^G^^  "  ^1^'  *""°'     *^>°  (83) 

(84) 
(85) 

(86) 
Using  Equation  (66)  the  Laplace  transform  will  be 

V  -  C^exp**^  +  Z^^Y.^"^^  (87) 


B.C. 8 

C  »  0, 

X  o  L, 

t>0 

B.C. 9 

C  »  0, 

0<X<L, 

t  « 

SC  T^v 

0<x<L, 

t>o 

where  n       t1/2 


=  [^J 


(88) 


Applying  the  boundary  conditions  to  eliminate  C-  and  C.  we  obtain 

k_L  8inh(L-x)q 
~   s(qLCo8hqL  +  k^LsinhqL)  (8^) 

b 

The  inverse  Laplace  transform  is 


J      /     kgLexp  8inh(L-x)qdB 


(90) 


7-ioo   8(qLCoshqL  +  k^LsinhqL) 

The  poles  are  s  »  0,  and  s  =  -^^  where  4  are  the  roots  of  the 
transcendental  equation 

Solution  of  Equation  (90)  by  residue  method  results  in 

kp(L-x) 
V  «-2 

«  j^2    (L-x)^n 

n^O  ^n^«M^  -^  ( V>^  V> 


r 

V  ,■ 


(  ) 
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-DC.dv 


Now   N  ,  -^   ^— 
x«»L     dx 


2 

kgPC^   ""  2D(kgL)\  Z,  '  -^^„  t 

1  +  k^L  "^     L 


y exp  '» 


.^  (93) 

Using  relationship  tan^  =  r-\  (94) 

n       kgL 

'=°'«n  =  <-»"    c.  ^^.  ..2,,y,                                     (95) 


I^TTTi^ITp 


-D^  ^t 
^n 


x=L       I  +  k^L  *  ^"•^^l  /  72     2 


H  

(96) 

k  DC  DC.     ,       k„       s 

L         G 


DC 
11m  N  a     i  .^^^ 

t-*    09  L  (98) 

Equation  (98)  Is  similar  to  Equation  (80). 

Equation  (96)  can  be  used  to  ascertain  the  time  required  to 

reach  steady-state  by  using  the  parameter  T  =  -^l  and  k  L,  thus  by 

-^  G 


Dt 
t*jr    uDjLu^}   I. lie  parameter  t  » 

the  previous  method 


n4i    CmV>'-'^g^ 


-D?  ^t 


(99) 

The  value  of  SST  for  different  values  of  T  and  k^L  are  shown  in 

G 

Figure  8. 
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CHAPTER  III 
EXPERIMENTAL  PROCEDURE 

A.  Experimental  Apparatus 
The  experimental  apparatus  consisted  of  a  diffusion  cell, 
pressure  measuring  system,  vacuum  system,  gas  cylinders  and  a  tempera- 
ture control  systeoi.  These  are  shown  In  Figure  9.  .. 

The  diffusion  cells,  shown  in  detail  In  Figures  10  and  11, 
consisted  of  a  stainless  steel  flange,  a  Pd-Ag  membrane,  and 
copper  gasket  held  together  by  four  bolts.  The  silver-palladium 
membranes,  1  1/8  inches  in  diameter,  were  punched  from  25%  silver- 
palladium  foil.  A  cavity  1/2  Inches  in  diameter  and  1/16  inches 
deep  was  machined  in  the  upstream  flange  to  serve  as  a  gas  chamber; 
a  gas  inlet  and  outlet  to  this  chamber  were  also  provided  as  shown. 
The  downstream  flange  was  similar  except  that  only  a  gas  outlet  was 
provided;  however,  a  thermocouple  well  was  machined  in  the  downstream 
flange.  The  downstream  cavity  was  connected  to  a  stainless  steel  (5641 cm^ 
tank  which  provided  enough  capacity  to  prevent  the  downstream  pressure 
from  changing  too  rapidly  during  measurement. 

A  thermocouple  gauge  was  used  to  measure  the  vacuum  on  the 
downstream  side  from  1-1000  microns.  A  Mcleod  gauge  was  used  to 
measure  the  pressure  change  in  the  downstream  side  of  the  apparatus 
during  the  diffusion  experiment.  A  pressure  regulator  was  used  to 
provide  the  desired  pressure  from  the  feed  gas  cylinder,  an  extra 
pressure  gauge  is  also  included  in  the  line. 

The  temperature  control  system  consists  of  a  temperature 
recorder-controller  and  a  tube  furnace.  A  chromel-alumel  thermocouple 
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was  used  to  measure  the  temperature  of  the  membrane. 

The  flow  through  the  system  was  controlled  by  a  toggle  valve, 
and  a  flowmeter  was  incorporated  In  the  upstream  side  to  regulate  the 
bleed  rate.  A  rupture  diaphragm  with  a  bursting  pressure  of  50  psig 
was  Installed  on  the  upstream  side. 

Two  Welch  vacuum  pumps  were  used  for  evacuating  the  system. 
Small  gas  cylinders  were  used  for  storage  of  feed  gas  and  for  the 
diffused  hydrogen  gas.  The  temperature  recorder  has  an  estimated  error 
of  ±1  C  and  the  Mcleod  gauge  has  an  estimated  error  of  ±0.05  mm  in 
the  range  of  pressure  used. 

B.  Experimental  Procedure 

The  transport  rate  of  hydrogen  through  the  palladium- silver 
membrane  was  measured  by  trapping  the  diffused  hydrogen  gas  in  the 
downstream  side  for  a  short  length  of  time  and  determining  the  pressure 
change  during  this  time  interval. 

The  experimental  procedure  is  important  because  a  different 
procedure  will  give  different  results.   In  all  these  experiments, 
unless  otherwise  specified,  the  first  set  of  experiments  is  always 
disregarded,  because  the  palladium- silver  membrane  gave  a-  lower 
transport  rate  when  used  for  the  first  time.   The  transport  rate  was 
found  to  increase  after  one  cycle  and  it  also  became  more  reproducible. 

Two  experimental  procedures  were  used  to  determine  the  trans- 
port rate:  Procedure  A  altered  the  pressure  while  keeping  the  tempera- 
ture constant.  In  Procedure  B  the  temperature  was  altered  while 
keeping  the  pressure  constant.   In  both  procedures,  the  pressure  and 
temperature  was  always  increased  from  a  lower  value  to  a  higher  one. 
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Altering  this  sequence  was  found  to  produce  different  results.  A 
more  detailed  description  of  these  procedures  is  given  in 
Appendix  B. 
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CHAPTER  IV 

EXPERIMENTAL  RESULTS  AND  CALCULATIONS 

A.  Experimental  Evaluation  of  the  Transport  Rp^«, 
Per  Unit  Area  

^  ^^  transport  rate  of  hydrogen  per  unit  area  through  0.106  nun 

25%  palladium-silver  foil  were  calculated  from  the  pressure  change  and 

time  data.  The  experimental  data  from  these  runs  are  tabulated  in 

Tables  1-19.  The  membrane  area  for  all  these  runs  was  1.268  cm^. 

Hydrogen  is  admitted  to  the  upstream  side  and  the  downstream 

pressure  increases  from  P^  to  P^  at  an  interval  time  0.  The  moles 

of  hydrogen  diffusing,  assuming  ideal  gas  behavior,  will  be: 

V!f-^) 


Cj 


''  ) 


"2  "  "l  "  RT   ^  (100) 

where   n^  =■  gm.  moles  at  the  start  of  the  experiment 
ng  «  gm.  moles  at  time  0 

P^  =  downstream  pressure  at  start  of  the  experiment  in  mm  Hg 
Pj  -  downstream  pressure  at  time  0  in  mm  Hg 
V^  "  volume  of  the  downstream  compartment  in  Cm^ 
T  ■  absolute  temperature  **K 
R  =»  gas  constant 
The  transport  rate  per  unit  area  will  be  given  by: 

n«  -  n    V  AP 

N.  ■  -^ i  a  _£ 

A     A0      A0RT  (101) 

where   N^  -  transport  rate  per  unit  area  in  gm.moles/cm^sec. 

V^  -  volume  of  downstream  compartment,  5641  cm^ 

A  =  area,  1.268  cm"' 

0  ■»  time,  sees. 
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^  »  Pj  -  Pjl^,  mm  Hg  . 

Substituting  the  corresponding  numerical  values,  Equation  (101) 
becomes: 

,  .  .  M  -, 5641 AP   ,  , ,     -2  AP 

L  ^  A   (1.268)(6.236  x  10^)  0T  "  ^'^^   '^  ^°   01  (^°2) 

Equation  (102)  is  the  working  equation  for  calculating  the 
transport  rate  per  unit  area  from  the  experimental  data.  The  values  of 
N^  are  also  tabulated  in  Tables  1-19. 

.  B.   Determination  of  the  Diffusivitv 

of  Hydrogen  in  Palladium-Silver  Alloy 

The  diffusivity  of  hydrogen  in  palladium- silver  alloy  may  be 

determined  from  Equation  (30)  and  the  diffusion  flux  from  Table  1  and 

t  Table  2. 

The  assumptions  for  these  calculations  are  (1)  Equation  (30) 
is  accurate  enough  for  the  purpose  of  these  calculations,  (2)  the 
adsorption  and  dissociation  process  is  fast  compared  to  diffusion  in 
the  solid  so  that  equilibrium  exists  between  the  gas  phase  and  the 
palladium- silver  surface. 

Using  Equation  (80)  therefore 


^A       T  (103) 


but 


^  „               gm. moles  H»       __       ^  „,   . 
C  »  r0  SiLiM 3- 2       gm.   at.Pd-A^ 

at.Pd-Ag         2  at.H  "^  107  gm.Pd-Ag 


jj  11.62  em.Pd-Ap; 
cm  Pd-Ag 


also  L  »  0.0102  cm 

r  -  0.36 


c 


78 


therefore 

^  (0.0102)(2)(107)   \  . 

(11.62)(0.36)  ^    e  ^^"^^ 

Equation  (104)  is  the  working  equation  for  determining  the 

dlffuaivity. 

From  the  tabulated  values  of  N.  in  Table  1  and  Table  2,  and 

A 

calculated  values  of  9   from  Equation  (30),  the  values  of  D  for  different 
temperatures  are  calculated  and  tabulated  in  Table  20  and  Table  21.  An 
average  value  of  D  is  taken  for  each  temperature  in  the  pressure  range. 
The  logarithm  of  these  values  are  plotted  against  1000/T  in  Figure  12. 
Figure  12  indicates  a  relationship  of  the  form: 

D  =  D^exp(-  -^)  (105) 

where   D  =  constant 
o 

AE^  =■  energy  of  activation  of  diffusion 
From  Figure  12  the  value  of  D  at  a  particular  temperature  is  found  to 
be  given  by: 

D  =  5.9  X  10"\xp(-  •^^)  (106) 

This  is  compared  with  an  accepted  value  (23)  of 

D  =-  4.3  X  10"^exp(-  ^^)  (107) 

for  pure  palladium  at  high  temperature;  although  D  is  roughly  the  same 
in  both  cases,  the  activation  for  diffusion  differs  in  the  two  materials. 

C.   Temperature  and  Pressure  Dependency  of  Flux 
of  Pure  Hydrogen  in  Palladium-Silver  Alloy 

The  fluxes  recorded  in  Tables  1  and  2  are  plotted  as  logN. 

versus  logP   .   A  straight  line  relationship  Is  found  as  shown  in 
mm 
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Figure  12.   Diffusivity  of  Hydrogen 
In  Palladium  Silver 
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Figure  13  and  Figure  15;  a  plot  of  logN^  versus  1000/T  also  gives  a 
straight  line  relationship  as  shown  In  Figure  14  and  Figure  16.  These  • 
results  suggest  that  the  flux  is  related  to  the  pressure  and  temperature 

by: 

1 

\   -  k^P"  exp  (.^  )  •         (108) 

where   k.  =■  constant 

AE  =  energy  of  activation 

—  »  constant 
n 

The  theoretical  verification  of  Equation  (108)  can  be  traced 
back  to  Equation  (80)  in  which 


DC 


but 


D  =  V^("Rr)  (110) 

C  =-k3P"exp(^-^  J  (HI) 


where   k-  =>  constant 

AHg  "  heat  of  solution 
The  equation  for  C  should  come  from  Equation  (30),  but  this 
can  be  approximated  by  Equation  (111).  Figures  17-18  show  that 
V_'  Equation  (111)  is  a  good  approximation  for  the  ranges  of  the  tempera- 

ture and  pressure  used  in  the  experiment.  Equation  (109)  therefore 
becomes: 

koD   r      •  AE„  -  AH, 


^A  "  T^  ^"^''P  -  (  RT J  "  kl"exp(-AE/RT)    (112) 

A  multiple  regression  analysis  technique  was  used  to  correlate 
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B:   3.79  X  10   gtn.  moles/cm  sec. 
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Figure  13.   The  Rate  of  Transport  of  Hydrogen 
Through  Palladium.  Silver  Membrane 
at  Different  Temperatures 
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B:   3.79  X  10   gm. moles/cm  sec, 

Feed  Gas  Composition:   Pure  H- 

Frocedure:   B 
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Figure  15.  The  Rate  of  Transport  of  Hydrogen 
Through  Palladium  Silver  Membrane 
at  Different  Temperatures 
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Figure  16.  Determination  of  the  Energy 
of  Activation 
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Figure  17.  A  Plot  of  Log  C  Versus  Log  P 
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Figure  18.  A  Plot  of  Log  C  Versus 
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Equation  (108)  with  the  data.   For  the  data  in  Table  1  and  Table  2 
this  resulted  in 

N^  =  2.,7.10-V-"e«p(^ij|^)  (114) 

when  the  feed  gas  is  pure  hydrogen. 

D.  Effect  of  Impurities  on  the  Transport 
of  Hydrogen  Through  Palladium!  Silver 

Table  3  to  Table  19  give  the  transport  rate  of  hydrogen 
through  palladium- silver  membranes  for  different  feed  gas  compositions 
and  various  bleed  rates.  Using  the  same  kind  of  plot  as  in  Figure  12 
and  Figure  13  of  the  previous  section  and  applying  similar  principles 
to  the  data  of  Tables  3  to  19,  resulted  in  the  data  shown  in  Figures 
19  to  52.   Except  for  carbon  monoxide  and  for  very  low  bleed  rate, 
the  experimental  data  for  hydrogen  diffusion  in  all  impurities  fitted 
Equation  (108).   The  results  of  the  least  square  fit  are  tabulated  in 
Table  22.  The  k^^  in  these  tables  indicate  that  total  pressure  was 
used  in  the  calculation,  and  k^*  when  partial  pressure  of  hydrogen  was 
used. 

From  the  values  indicated  in  Tables  3  to  16,  it  can  be  concluded 
that  bleeding  increases  the  transport  rate.   Equation  (108)  is  not 

suitable  for  expressing  the  relationship  between  transport  rate  and 

2 

bleed  rate.   It  was  found  that  a  plot  of  bleed  rate  B,  gm. moles/cm  sec, 

versus  B/N.  gave  a  straight  line  as  shown  in  Figures  53  to  77.  The 
relationship  is  given  by  the  equation: 
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B:   0.55  X  10   gm. mo lea /cm  sec, 


Feed  Gas  Composition:  30%  CH.,  10%  N,,  60%  H, 
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Figure  19.  Effect  of  Impurities  on  the  Rate 
of  Transport  of  lU  Through  Pd-Ag  Membrane 
at  Different  Temperatures 


fl 


0.9 


0.8 


^o    0-7 


00 

o 


0.6 


0.5 


0.4 


92 


1000 


2310 


mm 


2000  nun 
1690  mm 
1380  mm 


1069  mm 


1-3  1.4  1.5  1.6  1.7  1.8 


Figure  20.   Determination  of  the  Energy 
of  Activation 
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B:   1.08  X  10   gm. moles/cm  sec. 
Feed  Gas  Composition:   30%  CH,  ,  107.  N  ,  60%  H_ 
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Figure  21.  Effect  of  Impurities  on  the  Rate 
of  Transport  of  Hydrogen  Through 
Palladium  Silver  Membrane 
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Figure  22.   Determination  of  the  Energy 
of  Activation 
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Bt      1.88  X  10       gm. moles/ cm  sec. 


Peed  Gas  Composition:      30%  CH^,    10%  Nj,  60%  H, 
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Figure  23.  Effect  of  Impurities  on  the  Rate 
of  Transport  of  H-  Through  Pd-Ag  Membrane 
at  Different  Temperatures 
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Figure  24.     Determination  of  the  Energy 
of  Activation 


97 


-5  2 

B:  2.80  X  10   gm.moles/cm  sec. 


Feed  Gas  Composition:  30%  CH, ,  107.  N-,  607.  H 
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Figure  25.  E££ect  o£  Impurities  on  the  Rate 
of  Transport  of  Hg  Through  Pd-Ag  Membrane 
at  Different  Temperatures 
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Figure  26.  Determination  of  the  Energy 
of  Activation 
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B:  0.38  X  10   gm. moles /cm  sec. 

Feed  Gas  Composition:   10%  CH,  ,  107.  N  ,  80%  H^ 
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Figure  27.  Effect  of  Impurities  on  the  Rate 
of  Transport  of  H^  Through  Pd-Ag  Membrane 
at  Different  Temperatures 
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Figure  28.  Determination  of  the  Energy 
of  Activation 
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B:   0.52  X  10   gm. moles/cm  sec. 

Feed  Gas  Composition:  10%  CH. ,  107.  N.,  80%  H. 
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Figure  29.  Effect  of  Impurities  on  the  Rate 
of  Transport  of  H.  Through  Pd-Ag  Membrane 
at  Different  Temperatures 
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Figure  30.   Determination  of  the  Energy 
of  Activation 
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B:   1.11  X  10   gm. moles/cm  sec. 
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Figure  31.  Effect  of  Impurities  on  the  Rate 
of  Transport  of  H,  Through  Pd-Ag  Membrane 
at  Different  Temperatures 
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Figure  32.   Determination  of  the  Energy 
of  Activation 


105 


-5  2 
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Figure  33.  Effect  of  Impurities  on  the  Rate 
of  Transport  of  H-  Through  Pd-Ag  Membrane 
at  Different  Temperatures 


o 


1.1  r 


1.0 


0.9 


^o 


s 


o  0.8 


0.7 


0.6 


1.3 


106 


2 
3 
4 
5 


2000  mm 
1690  mm 
1380  mm 
1069  mm 


1.4 


1.5 


1.6 


1000 

T 


1.7 


1.8 


Figure  34.      Determination  o£  the  Energy 
of  Activation 
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B:   0.93  X  lO"  gm.mole/cm  sec. 
Feed  Gas  Composition:   25.67.  Ng,  74.47.  H, 


1.1  r 


C) 


1.0 


VO 
O 


O 


0.9 


0.8 


460°C 
420''C 
380°C 


300°C 


0.7 


c 


0.6 


3.0 


3.1 


3.2 
log  P 


nsn 


3.3 


3.4 


Figure  35.  Effect  of  Impurities  on  the  Rate 
of  Transport  of  H.  Through  Pd-Ag  Membrane 
at  Different  Temperatures 
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B:   2.46  X  10   gm. moles/ cm  sec. 

Peed  Gas  Composition:  25.6%  N2,  74.47,  H 
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Figure  37.  Effect  of  Impurities  on  the  Rate 
of  Transport  of  H^  Through  Pd-Ag  Membrane 
at  Different  Temperatures 
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Figure  38.  Determination  of  the  Energy 
of  Activation 
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B:   1.25  X  10   gin.  moles /cm  sec. 


Feed  Gas  Composition:   19.87.  CH,,  80.27.  H, 
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Figure  39.   Effect  of  Impurities  on  the  Rate 
of  Transport  of  H-  Through  Pd-Ag  Membrane 
at  Different  Temperatures 
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Figure  40.  Determination  of  the  Energy 
of  Activation 
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Figure  41.  Effect  of  Impurities  on  the  Rate 
of  Transport  of  H,  Through  Pd-Ag  Membrane 
at  Different  Temperatures 
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Figure  42.   Determination  of  the  Energy 
of  Activation 
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B:   0.86  X  10   gm. moles/cm  sec. 

Feed  Gas  Composition:   21.5%  CO.,  78.5%  H2 


c 


1.1 


1.0 


^o  °-9 


CO 

o 


0.8 


0.7 


0.6 


3.0 


3.1 


3.2 


3.3 


460°C 
420°C 

380°C 


340°C 


300°C 


3.4 


log  P 


mm 


Figure  43.  Effect  of  Impurities  on  the  Rate 
of  Transport  of  H_  Through  Pd-Ag  Membrane 
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Figure  44.  Determination  of  the  Energy 
of  Activation 
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B:   2.32  X  10   gm. moles /cm  sec. 
Feed  Gas  Composition:   21.5%  CO2,  78.57,  H^ 
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Figure  45.  Effect  of  Impurities  on  the. Rate 
of  Transport  of  H2  Through  Pd-Ag  Membrane 
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Figure  46.   Determination  o£  the  Energy 
of  Activation 


i»; 


119 

-5  2 

B:   0.99  X  10   gm. moles/cm  sec. 
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Figure  47.   Effect  of  Impurities  on  the  Rate 
of  Transport  of  H_  Through  Pd-Ag  Membrane 
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Figure  48.   Determination  of  the  Energy 
of  Activation 
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B:   2.52  X  10   gm. moles/cm  sec. 

Feed  Gaa  Composition:   21%  CO,  79%  H- 
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Figure  49.  Effect  of  Impurities  on  the  Rate 
of  Transport  of  H_  Through  Pd-Ag  Membrane 
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Figure  50.   Determination  of  the  Energy 
of  Activation 
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Figure  51.  Effect  of  Impurities  on  the  Rate 
of  Transport  of  H-  Through  Pd-Ag  Membrane 

at  Different  Temperatures 
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Figure  52.   Determination  of  the  Energy 
of  Activation 
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Figure  53.   A  Plot  of  B  Versus  B/N, 
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Figure  55.     A  Plot  of  B  Versus   B/N, 
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Figure  56.   A  Plot  of  B  Versus  B/N, 
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Figure  57.   A  Plot  of  B  Versus  B/N 
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10%  CH^,    107.  N2,   807.  H2 


r 


o 


3r 


o 

Q) 
U 


CO 
41 


o 


PQ 


Temperature: 

12  3  4 


300°C 


1 

s: 

2310 

nnn 

2 

a 

2000 

mm 

3 

s 

1690 

itmi 

4 

s 

1380 

mm 

5 

C3 

1069 

mm 

4 


b/n' 


Figure  58.  A  Plot  of  B  Versus  B/N, 
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Figure  59.   A  Plot  of  B  Versus  B/N 
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Figure  60.  A  Plot  of  B  Versus  B/N 
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Figure  62.  A  Plot  of  B  Versus  B/N. 


i 


136 


Feed  Gas « Compos  it  Ion:   25.6%  N,,  74. 47.  H, 
Temperature:   300°C 
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Figure  63.   A  Plot  of  B  Versus  B/N 
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Feed  Gas  Composition:   25.6%  H.,    74.4%  H2 
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Figure  64.   A  Plot  of  B  Versus  B/N. 
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Feed  Gas  Composition:      25.67,  Nj,    74.47.  H„ 
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Figure  65.     A  Plot  of  B  Versus   B/N 
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Figure  66.  A  Plot  of  B  Versus  B/N, 
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Feed  Gas  Composition:   25.6%  N2,  74.4%  Hj 
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Figure  67.  A  Plot  of  B  Versus  B/N, 
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Feed  Gas  Composition:   19.87.  CH,,  80.27.  H, 


Temperature:   300  C 
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Figure  68.  A  Plot  of  B  Versus  B/N, 
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Feed  Gas  Composition:   19.8%  CH,,  80.2%  H 
Temperature:   340°C 
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Figure  69.   A  Plot  of  B  Versus  B/N 
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Figure  70.  A  Plot  of  B  Versus  B/N 
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Figure  71.  A  Plot  of  B  Versus  B/N 
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Feed  Gas  Compos It ion:   19.8%  CH^,  80.2%  H^ 
Temperature:   460  C 
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Figure  72.     A  Plot  of  B  Versus   B/N 
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Feed  Gas  Composition:   21.57.  CO,,  78.5%  H 
Temperature:   300*'c 
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Figure  73.   A  Plot  of  B  Versus  B/N, 
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Feed  Gas  Composition:   21.57.  CO.,  78.57.  H_ 
Temperature:   340  C 
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Figure  74.   A  Plot  of  B  Versus  B/N. 
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Feed  Gas  Composition:      21.5%  CO.,    78.5%  H„ 
Temperature:      380°C 
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Figure   75.      A  Plot  of   B  Versus   B/N, 


U9 


o 


n 


3  r* 


Peed  Gas  Composition:   21.5%  CO,,  78.5%  H. 
Temperature:   420°C 
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Figure  76.   A  Plot  of  B  Versus  B/N, 


150 


Feed  Gas  Composition:   21.5%  CO^,  78.57.  H- 
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Temperature:      460°C 
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Figure    77.      A  Plot  of  B  Versus   B/N, 
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^  "  ^Amax  f-  ■*•  ^^2  (115) 

A 

Equation  (115)  is  similar  to  that  used  by  DeRosset  (6)  for 

constant  pressure  and  temperature.   The  significance  of  N     is 

Amax 

that  it  is  the  maximum  transport  rate  obtainable,  i.e.,  when  the  bleed 
rate  is  infinite.   Values  of  N^^  for  different  temperatures  and 
pressure  are  tabulated  in  Tables  23  to  27.  k^   signifies  the  effective- 
ness of  bleeding.   A  large  k^  means  that  bleeding  is  not  very  effective, 
while  a  small  k-  means  that  bleeding  is  quite  effective.  Values  of 
k^  for  different  temperatures  are  tabulated  in  Table  28. 
Equation  (115)  can  also  be  written  in  the  form: 

A       B 

N.    °  B  -  k„  (116) 

Amax        2 

Equation  (116)  will  be  used  to  illustrate  the  significance  of  k, . 

At  300  C  for  a  feed  composition  of  30%  CH.  ,  10%  N-  and  607.  H.  the 

value  of  kg  was  -3  x  lO'  .  On  substituting  this  value  into  Equation 

(116)  one  obtains: 

^'a        B 
_ ^j7j 

Amax   B  +  3  X  10 

A  plot  of  Equation  (117)  is  given  in  Figure  78.  This  shows 

V  ^^^^   *-^  the  value  of  k2  is  small,  the  transport  rate  will  approach  the 

maximum  value  even  for  a  small  bleed  rate.  To  find  B/N.  when  N  /N 

A       A  Amax 


o 


is  around  0.8,  for  this  particular  case.  Figure  78  gives 

-5  2 

B  =  1.2  X  10   gm. moles/cm  sec.   Then  for  a  pressure  of  1069  mm  I 

using  Figure  53  gives  B/N^  -  3.5.   This  means  that  to  obtain  80% 


c 
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TABLE  23 

EFFECT  OF  IMPURITIES  ON  THE  MAXIMUM  RATE  OF  TRANSPORT  ' 
OF  HYEROGEN  THROUGH  PALLADIUM  SILVER  MEMBRANE 

Feed  Gas  Composition:   30%  CH. ,  10%  N,,  60%  H- 

6  2 

Temp.  °C        Pressure.  iMi  He        "Amaii  '  ^°  8°.°olea/cni  sec. 


C") 


300 

1069 

4.3 

300 

1380 

5.0 

300 

1690 

5.8 

300 

2000 

6.2 

300 

2310 

7.0 

340 

1069 

5.3 

340 

1380 

6.4 

340 

1690 

7.2 

340 

2000 

7.6 

340 

2310 

8.1 

380 

1069 

6.2 

380 

1380 

6.9 

380 

1690 

7.5 

380 

2000 

7.9 

380 

2310 

8.6 

420 

1069 

6.6 

420 

1380 

7.6 

420 

1690 

8.2 

420 

2000 

8.6 

420 

2310 

9.4 

460 

1069 

7.2 

460 

1380 

8.5 

460 

1690 

9.1 

460 

2000 

9.6 

460 

2310 

10.6 

r 
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TABLE  24 


EFFECT  OF  IMPURITIES  ON  THE  MAXIMUM  RATE  OF  TRANSPORT 
OF  HYDROGEN  THROUGH  PALLADIUM  SILVER  MEMBRANE 


Feed  Gas  Composition; 


10%  CH^,  10%  Ng,  80%  Hj 


Temp.  C 


N. 


6  2 

X  10  gm .moles/cm  sec. 


c 


o 


300 

1069 

5.8 

300 

1380 

6.7 

300 

1690 

7.5 

300 

2000 

8.2 

300 

2310 

9.2 

340 

1069 

7.3 

340 

1380 

8.0 

340 

1690 

8.9 

340 

2000 

9.8 

340 

2310 

10.6 

380 

1069 

7.7 

380 

1380 

8.7 

380 

1690 

9.7 

380 

2000 

10.5 

380 

2310 

11.4 

420 

1069 

8.2 

420 

1380 

9.3 

420 
420 

1690 
2000 

10.2 
10.9 

420 

2310 

12.0 

460 

1069 

8.4 

460 

1380 

9.8 

460 

1690 

10.4 

460 
460 

2000 
2310 

11.0 
12.5 
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TABLE  25 


r 


EFFECT  OF  IMPURITIES  ON  THE  MAXIMUM  RATE  OF  TRANSPORT 
OF  HYDROGEN  THROUGH  PALLADIUM  SILVER  MEMBRANE 


Feed  Gas  Composition: 


o 


TemPi  °C 

Pressure,  nm  Hr 

300 

1069 

300 

1380 

300 

1690 

300 

2000 

300 

2310 

340 

1069 

340 

1380 

340 

1690 

340 

2000 

340 

2310 

380 

1069 

380 

1380 

380 

1690 

380 

2000 

380 

2310 

420 

1069 

420 

1380 

420 

1690 

420 

2000 

420 

2310 

460 

1069 

460 

1380 

460 

1690 

460 

2000 

460 

2310 

25 

.6% 

N^,  74.4% 

«2 

N,    X 

Amax 

10« 

gm.mo 

les/cm  sec. 

5.4 
6.2 
6.9 
7.5 
8.6 

• 

6.4 
7.4 
7.9 
8.4 
9.0 

6.7 
8.0 
8.7 
9.4 
10.1 

7.6 

8.7 

9.3 

10.0 

10.9 

8.5 

9.1 

10.0 

10.9 

11.7 

JLJUti't 
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TABLE  26 


EFFECT  OF  IMPURITIES  ON  THE  MAXIMUM  RATE  OF  TRANSPORT 
OF  HYDROGEN  THROUGH  PALLADIUM  SILVER  MEMBRANE 


C" 


o 


F 

eed  Gas 

Composition:   19.8% 

CH^,  80.2%  H2 

• 

Temp.  °C 

Pressure,  mm  Hj^ 
1069 

N.    X  10^ 
Amax 

2 
gm. moles/cm  sec. 

300 

5.5 

300 

1380 

6.1 

300 

1690 

6.9 

300 

2000 

7.7 

300 

2310 

8.5 

3A0 

1069 

6.4 

340 

1380 

7.5 

340 

1690 

8.1 

340 

2000 

8.9 

340 

2310 

10.0 

380 

1069 

7.5 

380 

1380 

8.6 

380 

1690 

9.4 

380 

2000 

10.2 

380 

2310 

10.9 

420 

1069 

8.5 

420 

1380 

9.0 

420 

1690 

10.1 

420 

2000 

11.0 

420 

2310 

12.0 

460 

1069 

8.4 

460 

1380 

9.7 

460 

1690 

10.7 

460 

2000 

11.6 

460 

2310 

12.7 

r 
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TABLE  27 

EFFECT  OF  IMPURITIES  ON  THE  MAXIMUM  RATE  OF  TRANSPORT 

OF  HYKIOGEN  THROUGH  PALLADIUM  SILVER  MEMBRANE 

Feed  Gas  Composition:  21.57.  CO-,  78.5%  ^o   • 

6  2 

_     o_  __          „  N.    X  10  gm. moles/cm  sec. 

Temp.  C  Pressure t  mm  Kg  Amax [ 

300  1069  5.8 

300  1380  6.5 

300  1690  7.0 

300  2000  7.5 

300  2310  8,1 

340  1069  6.2 

340  1380  7.5 

340  1690  8.4 

340  2000  9.2 

r'--^                             340  2310  10.3 

380  1069  7.5 

380  1380  8.8 

380  1690  9.6 

380  2000  10.2 

380  2310  10.7 

420  1069  8.0 

420  1380  9.1 

420  1690  9 . 7 

420  2000  10.5 

420  2310  11.3 

460  1069  8.4 

460  1380  9.6 

460  1690  10.4 

460  2000  11.0 

460  2310  12.0 


c 
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o 


o 


TABLE  28 


TABULATION  OF  k^   IN  EQUATION  (115) 
FOR  DIFFERENT  GASES  AND  TEMPERATURE 


Feed  Gas  Composition 
307.  CH^,  10%  N^,  60%  H2 


10%  CH^,  10%  N2,  80%  H2 


25.6%  Ng,  74.4%  H2 


19.8%  CH^,  80.2%  H2 


21.5%  CO2,  78.5%  Hj 


Temperature  C 

-(k2  X  10^) 

300 

3.0 

340 

1.4 

380 

0.8 

420 

1.0 

460 

2.0 

300 

0.8 

340 

0.8 

380 

1.0 

420 

0.8 

460 

0.8 

300 

0.4 

340 

0.5 

380 

0.6 

420 

0.8 

460 

0.6 

300 

1.0 

340 

0.8 

380 

0.8 

420 

0.8 

460 

1.0 

300 

0.8 

340 

0.6 

380 

0.6 

420 

0.4 

460 

0.4 

i 


C 


i 


c 


1.0  r 
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5  2 

B  X  10     gm. moles/cm  sec. 


Figure  78.     A  Plot  of 


N. 


^Vmx 


Versus  B 


n 
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efficiency  by  bleeding,  the  bleed  rate  must  be  3.5  times  larger  than 
the  transport  rate,  i.e.,  almost  78%  of  the  feed  must  be  discharged  in 
the  bleed. 

The  theoretical  verification  of  Equation  (116)  can  be  traced 
to  Equation  (97): 


°°i  /   ^G   \ 


if 

DC 


L     Amax 
than 


-  -  Na„..  (119) 


B        ^G 

T^TT  ' 1  <i20) 


G   L 
From  Equation  (120)  therefore 


In  other  words  k_  is  directly  proportional  to  the  bleed  rate. 
The  significance  of  k  is  that  it  is  a  hypothetical  mass  transfer 
coefficient  in  the  gas  phase,  where  the  concentrations  in  the  gas  phase 
are  expressed  in  terms  of  the  concentration  of  the  solid. 

If  Equation  (119)  is  true,  this  means  that  the  effect  of 
(~)  impurities  la  purely  in  the  form  of  gas  film  resistance.   This  was  not 

found  to  be  true.  When  the  values  of  N.    are  fitted  to  Equation  (108), 
the  usual  N^^^,  T,  P  plot  of  which  are  shown  in  Figures  79  to  88,  one 
gets  the  fitted  constant  tabulated  in  Table  29.  The  energy  of  activation 
was  higher  when  impurities  wereintroduced,  which  indicates  that  there 
is  a  surface  effect  such  as  adsorption  which  is  causing  this  deviation.' 
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Feed  Gas  Composition:   307.  CH, ,  10%  N_,  607.  H 
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Figure  79.   Effect  of  Impurities  on  the  Maximum 

.  Rate  of  Transport  of  H2  Through  Pd-Ag  Membrane 

at  Different  Temperatures 
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Figure  80.   Determination  of  the  Energy 
of  Activation 
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Feed  Gas  Composition:   107.  CH. ,  10%  N2,  80%  H- 
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Figure  81.  Effect  of  Impurities  on  the  Maximum 

Rate  of  Transport  of  H2  Through  Pd-Ag  Membrane 

at  Different  Temperatures 
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Figure  82.   Determination  of  the  Energy 
of  Activation 
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Feed  Gas  Composition:   25.6%  Ng,  74.4%  H. 
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Figure  83.  Effect  of  Impurities  on  the  Maximum 
Rate  of  Transport  of  Hg  Through  Pd-Ag  Membrane 
at  Different  Temperatures 
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Figure  84.  Determination  of  the  Energy 
of  Activation 
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Feed  Gas  Composition:   19.8%  CH,,  80.2%  H„ 
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Figure  85.   Effect  of  Impurities  on  the  Maximum 
Rate  of  Transport  of  H2  Through  Pd-Ag  Membrane 
at  Different  Temperatures 
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Figure  86.   Determination  of  the  Energy 
of  Activation 
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Feed  Gas  Composition:   21.5%  CO2,  78.5%  H 
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Figure  8<7.  Effect  of  Impurities  on  the  Maximum 

Rate  of  Transport  of  H-  Through  Pd-Ag  Membrane 

at  Different  Temperatures 
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Figure  88.   Determination  of  the  Energy 
of  Activation 
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TABLE  29 


MULTIPLE  REGRESSION  OF  THE  CONSTANT  IN  EQUATION  (108) 
FOR  THE  MAXIMUM  TRANSPORT  RATE 


Feed  Gas  Compositior 

I 

\ 

^ 

1 
n 

0.50 

AE 
R 

307,  CH^,  107,  N2,  607. 

»2 

1.1 

X  10"^ 

8.5 

X  10"^ 

-1160 

107,  CH^,  107,  N^,  807, 

«2 

8.0 

X  10"^ 

7.2 

X  10'^ 

0.51 

-  850 

25.67.  N2,  74.47,  H^ 

1.1 

X  10"^ 

9.4 

X  lO"'' 

0.48 

-  990 

19.87,  CH^,  80.27,  H^ 

1.1 

X  10"^ 

9.8 

X  10"^ 

0.52 

-1140 

21.57,  CO2,  78.57,  H^ 

1.2 

X  10-^ 

1.1 

X  10"^ 

0.48 

-1000 
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The  observed  value  of  the  energy  of  activation  shows  that  carbon 
monoxide  has  the  largest  energy  of  activation  followed  by  methane, 
carbon  dioxide  and  nitrogen.  In  Chapter  II  it  was  pointed  out  that 
gases  which  have  the  largest  heat  of  adsdrption  have  the  longest 
residence  time  at  the  surface.  The  longer  the  residence  time  of  the 
impurities  in  the  surface,  the  greater  will  be  the  energy  of  activa- 
tion. Therefore  if  the  heats  of  adsorption  of  the  impurities  were 
known,  the  energy  of  activation  of  the  transport  process  would  fall 
in  the  same  sequence  as  the  heats  of  adsorption.   It  is  also  observed, 
as  might  be  expected,  grease  lowers  the  transport  rate. 

E.   Time  to  Reach  Steady  State 
The  time  to  reach  807.  steady  state  in  the  experiment  varies 
from  a  few  seconds  to  less  than  two  minutes .   This  is  in  the  range  that 
is  expected  by  using  Equations  (81)  and  (99). 


r 
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CHAPTER  V 

DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 

A.  Transport  Rate  o£  Pure  Hydrogen  Through 
Palladium- Silver  Membranes 

The  transport  rates  £ound  in  this  study  are  several  times 
greater  than  those  observed  by  most  earlier  investigators.  This  is 
to  be  expected  for  several  reasons.  In  Chapter  II  it  was  pointed 
out  that  in  the  Orphase  region  hydrogen  is  more  soluble  in  palladium- 
silver  alloys  than  in  pure  palladium.  The  reverse  is  true  for  the 
/3-phase  region.  Hurlbert  and  Konecny  (5)  employed  pure  palladium 
membranes  whereas  palladium- silver  membranes  were  used  in  this  inves- 
tigation. Since  all  of  this  work  was  done  under  conditions  where  the 
hydrogen  concentrations  were  relatively  small,  it  is  easily  seen  that 
the  hydrogen  transport;  rate  through  the  membrane  should  be  higher  for 
the  alloy  than  for  pure  palladium.  Kydonlus  (33)  did  not  bleed  his 
gas,  and  it  seems  likely  that  his  transport  rates  were  lowered  by  the 
presence  of  accumulated  impurities  on  the  upstream  side  of  the  membrane. 

The  experimental  transport  rate  data  correlate  very  well  with 

an  equation  containing  a  term  p     where  —  has  an  average  value  of 

0.62.  This  is  approximately  the  average  of  the  values  of  the  slopes  of 

the  isotherms  shown  in  Figure  17.   Kydonlus,  using  the  same  kind  of 

membrane  but  at  lower  temperatures,  obtained  a  value  of  —  about  0.70. 

n 

The  energy  of  activation  found  In  this  work  Is  lower  than 
that  found  by  previous  Investigators.   Their  values  range  from  2.5  to 
36  kcal/mole.  The  energy  of  activation  based  on  the  transport  rate  is 
an  overall  activation  energy,  and,  if  adsorption  or  dissociation 
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contributes  to  the  transport  process  (such  as  is  the  case  when  impurities 
are  present)  AE  should  be  greater  than  that  for  diffusion  alone.   Thus, 
large  energies  of  activation  strongly  suggest  that  surface  processes 
are  rate  limiting. 

Bleeding,  therefore,  appears  to  accomplish  two  things:  firstly, 
it  dilutes  the  impurities  present  at  the  upstream  face  of  the  membrane; 
secondly,  it  should  result  in  a  lower  energy  of  activation  since  it 
would  tend  to  reduce  the  number  of  adsorption  sites  occupied  by 
impurity  molecules.   There  may  be  a  thermal  effect  also.   The  incoming 
feed  gas  is  approximately  at  room  temperature  while  the  palladium- 
silver  membrane  is  at  a  higher  temperature,  and  the  difference  in 
temperature  would  be  sufficient  to  cause  convection.   Moreover,  since 
the  cooler  molecules  have  lower  velocities,  the  rate  at  which  they 
strike  unit  area  of  surface  would  also  be  decreased  as  compared  to  a 
hot  gas.   However,  the  diffusion  rate  in  the  gas  phase  is  so  high  that 
this  effect  appears  to  be  negligible. 

The  calculation  of  the  diffusivity  for  the  palladium- silver 
gives  a  value  which  is  slightly  lower  than  that  for  pure  palladium. 
This  is  to  be  expected  since  hydrogen  does  not  diffuse  rapidly  through 
silver. 

B.   Effect  of  Impurities  on  the  Transport  Rate 
of  Hydrogen  Through  Palladium-Silver  Membrane 

As  shown  in  Figures  19  to  52,  the  introduction  of  impurities 

into  hydrogen  lowers  the  transport  rate.   This  can  be  attributed  to 

the  following  effects:   (1)  hydrogen  is  diluted  by  the  impurities; 

(2)  impurities  hinder  the  diffusion  of  hydrogen  in  the  gas  phase; 
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(3)  occupation  of  adsorption  and  dissociation  sites  by  impurities 
result  in  a  reduction  of  hydrogen  adsorption  and  dissociation.   The 
overall  effect  of  impurities  is,  therefore,  to  lower  the  concentration 
of  hydrogen  atoms  in  the  upstream  surface  of  the  palladium- silver 
membrane,  and  since  the  diffusion  rate  in  the  alloy  is  directly  pro- 
portional to  concentration,  the  transport  rate  decreases. 

It  is  natural  to  assume  that  good  mixing  occurs  in  the  upstream 
cavity  of  the  diffusion  cell.   If  this  is  the  case,  and  there  are  no 
surface  effects,  then  Equation  (C-8)  of  Appendix  C  should  hold.   This 
was  found  not  to  be  true,  since  Equation  (115)  rather  than  Equation 
(C-8)  fits  the  data.  This  implied  that  two  regions  exist  in  the 
/"  ^  upstream  side,  (1)  a  region  where  plug  flow  exists,  and  (2)  a  region 

near  the  membrane  where  an  impurity-rich  layer  is  formed.  This  layer 
constitutes  a  resistance  to  diffusion  from  the  gas  phase  to  the 
surface  of  the  membrane  and  leads  to  the  postulation  of  a  hypothetical 
mass  transfer  coefficient,  k^,  based  on  gas  film  resistance.   The  mass 
/  v^ transfer  coefficient  was  found  to  be  directly  proportional  to  the  bleed 
rate,  as  shown  by  agreement  of  the  data  with  Enuation  (121).   Since 
the  mass  transfer  coefficient  is  inversely  proportional  to  the  gas 
film  resistance,  bleeding  the  feed  gas  reduces  the  gas  film  resistance. 

The  multiple  regression  analysis  of  Equation  (108)  for  N.  and 
^Amax  ^^°"^  ^^^^   n  ^^  generally  smaller  when  impurities  are  introduced; 
an  exception  is  grease.  This  is  possibly  caused  by  the  occupation  of 
adsorption  sites  for  hydrogen  by  impurities.  Gases  such  as  CH^,N_,  and 
COg  appear  to  be  physically  adsorbed  on  palladium-silver;  gases  such  as 
CO  appear  to  be  chemisorbed,  and  may  occupy  most  of  the  dissociation 
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sites  for  hydrogen.   It  can  be  deduced  from  Chapter  II  that  these 
gases  have  a  larger  residence  time  on  the  membrane  than,  hydrogen.  Since 
Increasing  the  total  pressure  increases  the  concentration  of  both 
hydrogen  and  impurities  in  the  gas  phase,  and  since  the  impurities  have 
a  longer  residence  time  than  hydrogen,  the  effect  of  the  impurities 
is  greater  than  one  would  predict  from  low  mole  fraction,  and  greater 
at  high  pressures  than  at  low  pressures.  This  lowers  the  value  of 

—  when  compared  with  a  feed  gas  of  pure  hydrogen. 

The  multiple  regression  analysis  for  N.  and  N.    also  shows 

A      Amax 

that  the  energy  of  activation  for  diffusion  is  higher  when  impurities 
are  introduced.   Since  the  energy  of  activation  is  based  on  the  overall 
transport  rate,  it  is  really  an  overall  activation  energy.   If  the 
impurities  adsorb  at  surface  of  the  membrane,  it  would  increase  the 
energy  of  activation  above  that  observed  for  pure  hydrogen,  and  it 
is  to  be  expected  that  a  gas  which  chemlsorbs  will  have  a  higher 
energy  of  activation  than  one  which  is  physically  adsorbed.  Gases 
with  unsaturated  double  bonds  were  found  to  have  a  higher  energy  of 
activation  for  diffusion  than  chemically  inactive  gases,  such  as  N_. 
The  transport  rate  of  hydrogen  containing  carbon  monoxide  as 
an  impurity  is  greatly  reduced  when  the  temperature  is  between  300°C 
to  380  C;  at  higher  temperatures,  i.e.,  from  420  to  460°C  this  system 
behaves  essentially  the  same  as  systems  with  inactive  gas  impurities. 

The  large  decrease  of  the  transport  rate  in  the  temperature  range 

o       ° 
300  C  to  380  C  suggests  a  poisoning  effect  of  the  carbon  monoxide 

similar  to  that  of  hydrogen  sulfide  observed  by  other  investigators  (5). 
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It  also  indicates  that  most  of  the  adsorption  and  dissociation  sites 
for  hydrogen  are  occupied  by  carbon  monoxide.   However,  sulfide  poison- 
ing is  stronger,  since  it  permanently  damaged  the  membrane  once 
^^  poisoned  by  hydrogen  sulfide.   Whereas,  the  poisoning  effect  of  carbon 

monoxide  is  removed  by  raising  the  temperature  to  420°C.   This  suggests 
that  the  bond  formed  between  palladium  and  carbon  monoxide  at  the 
low  temperatures  is  unstable  at  higher  temperatures. 

C.   Conclusion 

The  equation  for  the  transport  of  hydrogen  through  palladium, 
i.e..  Equation  (112),  fits  the  experimental  data  quite  well,  but 
the  transport  rates  observed  differ  by  a  factor  of  2  or  3  from  those 
(      )  obtained  by  recent  investigators.   Compared  with  some  investigators 

in  the  past,  for  which  the  transport  rate  varied  as  much  as  1000  times, 
this  agreement  is  fairly  reasonable. 

Equation  (112)  is  an  approximation  of  a  more  complicated 
equation,  but  is  reasonably  good  for  a  given  range  of  conditions. 

Introduction  of  Impurities  lowers  the  transport  rate  either  by 
the  dilution  effect  of  the  nondif fusing  component  in  the  gas  phase  or 
by  a  surface  effect,  or  both.   The  surface  effect  is  caused  by  the 
reduction  of  adsorption  sites  and  dissociation  sites  of  hydrogen  by 
the  impurities,  resulting  in  a  reduction  of  the  absorbed  hydrogen 
atom  concentration  just  Inside  the  upstream  membrane  surface.   If  the 
Surface  effect  is  so  large  that  it  becomes  the  rate- limiting  process 
in  the  transport  of  hydrogen  through  palladium,  it  is  said  that  the 
impurities  poison  the  membrane.   This  is  probably  caused  by  the 
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chemlBorptlon  of  impurities  at  the  surface  or  by  the  chemical  reaction 

of  the  impurities  with  the  surface  atoms.  Carbon  monoxide  and  grease 

poison  the  membrane  over  the  temperature  range  of  300°C  to  380**C,  but 

appear  not  to  do  so  at  temperatures  above  A20°C.  If  the  surface 

effect  is  small,  such  as  in  the  case  of  CH.,  CO.,  and  N.,  it  is 

believed  that  only  physical  adsorption  of  the  Impurities  occur  at  the 

surface. 

When  impurities  are  introduced  and  the  feed  gas  is  bled, 

Equation  (115)  fits  the  experimental  data  quite  well.  This  equation 

combined  with  Equation  (112)  for  expressing  N     fit  the  experimental 

data  for  a  wide  range  of  bleed  rate,  temperature,  and  pressure. 

Equation  (115)  suggested  that  there  exists  an  Impurity  rich  gas  film  near 

the  surface  of  the  membrane.   Comparing  H,    with  that  for  the  trans- 

Amax 

port  rates  of  pure  hydrogen  on  the  same  basis  and  conditions  indicates 

that  N^^^^  is  lower  than  that  for  pure  hydrogen.  Furthermore,  the 

multiple  regression  analysis  of  N.    shows  that  the  value  of  —  is 

Amax  n 

lower  and  the  energy  of  activation  is  higher  than  that  for  pure  hydrogen. 
All  this  suggests  that  surface  effect  is  present  when  impurities  are 
introduced . 


CHAPTER  VI 

RECOMMENDATIONS  FOR  FUTURE  STUDIES 

In  this  study,  no  positive  confirmation  as  to  the  exact  nature 
mji  of  the  reaction  between  the  metal  and  poisoning  gaf,  has  been  established. 

In  the  light  of  the  development  of  recent  analytical  methods,  it  may 
be  possible  to  determine  the  nature  of  the  bond  between  the  metal  and 
poison. 

The  effect  of  silver  content  on  the  diffusivity  can  also  be 
examined  more  extensively.   This  would  lead  possibly  to  the  study  of 
the  order-disorder  phenomena  in  the  palladium-silver  alloy.   The 
assumption  in  this  work  is  that  palladium-silver  atoms  are  completely 
<-"  .  random  or  disordered.   An  investigation  of  the  effect  of  other  impuri- 

ties  than  those  studied  here  is  also  recommended. 

The  effect  of  cell  construction,  anealing  procedures  of  the 
membrane,  and  repeated  cycling  to  determine  membrane  stability  and 
transport  efficiency  are  also  some  of  the  areas  which  need  further 
study. 
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NOMENCLATURE 

a  lattice  conBtant;  constant  In  Equation  (27) 

a  (T)  the  partition  function  for  the  internal  degree  of 
freedom  of  an  absorbed  hydrogen  atom  referred  to 
the  ground  state  at  infinite  separation 

A  area  of  membrane 

B  bleed  rate 

C.C. ,C  concentration 
1  o 

C.  concentration  of  A 
A 

^l'^2  constant  of  integration  in  Equation  (74) 

D  dlffusivity 

D  constant  defined  by  Equation  (53) 

E„  e„  in  the  molal  basis 

^uu  ^ou  i"  *^h6  molal  basis 
(In      nn 

E_  e_  in  the  molal  basis 

AE  energy  of  activation 

AEjj  energy  of  activation  of  diffusion 

E(X-Y)  bond  energy  between  X  and  Y 

f  geometric  factor 

F  Helmholtz  free  energy 
,f 


^)  AF     free  energy  of  formation 


AF     free  energy  of  activation 

g(Njj,Nj^)    total  number  of  distinguishable  configuration  of  N 
atoms  " 

h    Planck  constant 

AH    heats  of  adsorption 

AH_    heats  of  solution 
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k  Boltzmann  constant 

k. ,k. jk-.k-.k.  constants 

k^*,k  hypothetical  mass  transfer  coefficient 

K  equilibrium  constant 

K(T),K'(T)    constant  defined  by  Equation  (17);  K(T)  in  molal 
basis 

L  thickness  of  membrane 

M  mean  molecular  weight  per  constituent  atom 

Mj,  atom  mass  of  hydrogen 

M  moment  of  inertia  of  the  hydrogen  atom 

n(x)    the  number  of  diffusing  particles  per  cm  on  the 
plane  located  at  x  at  the  instant  t 

Hj^.ng    gram  moles  of  hydrogen  at  the  start  and  at  time  0 
respectively 

N.    molal  flux  of  A 

N„    number  of  hydrogen  atoms  in  the  lattice 

N-    number  of  interstitial  sites 

Njjjj    number  of  pairs  of  neighboring  hydrogen  atoms  in 
a  particular  configuration 

N      defined  by  Equation  (8) 

as 

N    defined  by  Equation  (2) 

P,P    pressure 

Pf  2)    equilibrium  pressure  of  the  a-^  phase  transition 

^t'^f    downstream  pressure  at  the  start  of  the  experiment 
and  at  time  0  respectively 

AP    P^  -  P 

P^    probability  for  a  given  interstitial  hydrogen  to 
make  a  Jump  per  second 

q    defined  by  Equation  (75) 
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4jj(T)    partition  function  for  vibration  of  absorbed  H 
atom  relative  to  the  lattice 

q^    the  initial  heat  of  adsorption 

Q    partition  function 

+ 
Q     partition  function  of  the  normal  state 

Q  partition  function  of  the  activated  state 

R  gas  constant 

R^  chemical  reaction  rate  of  A 

AS  entropy  of  activation 

SST    fraction  of  the  steady  state  value  approach  at 
dimension less  time  T 

t  time 

T  °  temperature 

T  ,T  melting  temperature 

T  critical  temperature 

V  defined  by  Equation  (66) 

* 

V     molar  average  velocity 

V  molar  volume;  Laplace  transform  of  v 
X.     weight  fraction  of  silver 

z    number  of  nearest  neighbors 

Greek  Letters 

/3    constant 

e^jj    total  interaction  energy  at  saturation  per  atom 

'^H    ^^^   energy  required  to  remove  an  absorbed  H  atom, 
far  removed  from  any  other  absorbed  H  atom,  from 
its  lowest  state  in  the  metal  to  rest  at  infinite 
dispersion 
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e_    energy  required  to  dissociate  a  hydrogen  molecule 
from  its  lowest  state  Into  two  free  atoms  In  its 
lowest  state 

0         fraction  of  Interstial  sites  occupied 

X„    absolute  activity 
H 

X    distance  between  two  atomic  planes  where  diffusing 
is  occuring 

]U  chemical  potential 

jU  electronegativity 

V  frequency  of  vibration  of  linear  oscillator 

p  spin  weight  of  the  proton 

T  adsorption  residence  time 

T    dimensionless  time  in  Equation  (81)  and 
Equation  (99) 

T     constant  defined  by  Equation  (36) 

0    time 
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APPARATUS 
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Diffusion  Cell:   Type  316  Stainless  Steel 

Chemical  Engineering  Department,  University  of 
Florida,  Gainesville,  Florida 
Electric  Furnace:   1400  Watts.  Maximum  Temperature  1950  F 

Fisher  Scientific  Company,  Pittsburgh,  Pennsylvania 
Thermocouple  Vacuum  Gauge:   Type  GTC-100.  0-1000  Microns  Range 

Consolidated  Vacuum  Corporation 
Rochester ,  New  York 
Temperature  Controller:   CHR.  AL.  Thermocouple,  Series  60  Control 

Leeds  and  Northrup,  Philadelphia,  Pennsylvania 
Timer:   Range  0.1-9999.9  Seconds 

Precision  Scientific  Co.,  Chicago,  Illinois 
Vacuum  Gauge:   Range  0.0001-25  mm,  Mcleod  Vacuum  Gauge 

Todd  Scientific  Co.,  Springfield,  Pennsylvania 
Duo- Seal  Vacuum 

Pump:   Welch  Series  1400B,  Duo-seal  Vacuum  Attainable,  0.1  Micron 
W.  M.  Welch  Manufacturing  Company,  Chicago,  Illinois 
Safety  Relief  Valve:  Model  No.  21,  Max.  Pressure  50  pslg 

The  Matheson  Company,  Inc.,  Bast  Rutherford, 
New  Jersey 
Toggle  Valve:  Model  1254 

Crawford  Fitting  Co.,  Cleveland,  Ohio 
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Metering  Valve:  Model  2PY  280  Metering  Valve  with  Vernier  Handle 

Hoke,  Inc.,  Cressklll,  New  Jersey 
Flow  Meter:   Model  995  Bantam  Flo-Gages-5-100  cc/mm  of  Carbon  Dioxide 
h  Hoke,  Inc.,  Cressklll,  New  Jersey 

Pressure  Gauge:   30"-100  pslg 

U.  S.  Gauge,  Sellersvllle,  Pennsylvania 
Fittings  and  Tubing:   1/8"  to  1/4"  fittings  and  tubing 

Crawford  Fitting  Co.,  Cleveland,  Ohio 
Silver  Palladium  Foil:   6"  x  6"  x  0.004"  Pd.  25%  Ag  Foil 

Engelhard  Industries,  Newark,  New  Jersey 
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APPENDIX  B 

EXPERIMENTAL  PROCEDURE 

A.  Preparation  of  the  Palladium-Silver  Membrane 

1.  Soak  the  palladium-silver  membrane  and  copper  gaskets 
in  ether  for  one  hour,  rinse  twice  with  distilled  water  and  drip. 

2.  Mount  the  palladium- silver  membrane  and  copper  gaskets 
in  the  diffusion  cell. 

3.  Evacuate  both  sides  of  the  membrane  by  opening  valve 
V4,  V5,  V7,  V8,  V9,  and  Vll.  All  other  valves  being  closed.  Heat  to 
460°C  for  one  hour  while  being  continuously  evacuated.  Cool  to  300  G 
while  continuously  evacuating. 

4.  When  a  temperature  of  300**C  is  reached,  close  valve  V4  and 
open  valve  V2  with  pressure  set  at  1069  mm  Hg.  Then  adjust  valve  V6 

so  that  it  reads  a  specific  bleed  rate  in  the  flowmeter.  Maintain 
.  this  condition  for  one  hour. 

5.  Change  pressure  to  2310  mm  Hg  and  maintain  other 
variables  at  same  valves  for  one.  hour. 

6.  Change  pressure  to  1069  vm  Hg  and  heat  to  460  C  for  one 
hour.  All  other  variables  maintain  at  the  same  valves. 

(^    )  7.   Repeat  step  5. 

8.   The  palladium-silver  membrane  is  now  ready  for  the 
experiment.  .If  in  subsequent  use  the  bolts  become  loose,  the  cell  is 
dismounted  at  both  connectors  and  the  bolts  are  tightened. 

B.  Procedure  A 
1.   Evacuate  both  sides  of  the  membrane  and  heat  to  460  C  for 
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one  hour  while  continuously  evacuating. 

2.  Cool  to  300  C  while  continuously  evacuating. 

3.  When  the  temperature  reaches  300°C,  close  valve  V4,  open 
valve  V2  with  the  pressure  set  at  1069  nm  Hg,  and  adjust  valve  V6  so 
that  the  flowmeter  reads  the  desired  bleed  rate.  Maintain  this  condi- 
tion for  one  hour. 

4.  Open  valve  VIO  for  5  minutes  and  take  reading  of  the 
initial  pressure  by  the  Mcleod  Guage.   Set  timer  and  close  valve  VIO 
and  Vll.  After  a  suitable  time  interval,  close  valve  V7  and  open  valve 
VIO;  take  reading  of  the  final  pressure  and  time.   The  corrected  final 
pressure  is  then  obtained  from  Figure  B-  1  .   (This  is  necessary  since 
the  gas  has  expanded  to  the  volume  occupied  by  the  Mcleod  Gauge.) 

5.  Close  valve  VIO,  open  valve  V7  and  Vll,  and  change  to  the 
next  higher  pressure  and  adjust  to  the  desired  bleed  rate.   After 

15  minutes  repeat  step  4. 

6.  Repeat  step  5  until  the  pressure  reaches  2310  mm  Hg; 

then  change  pressure  to  1069  ram  Hg  and  heat  to  the  next  higher  tempera- 
ture for  one  hour. 

7.  Finally,  repeat  step  4  to  step  6  until  temperature  of 
460  C  is  reached. 

C.   Procedure  B 

1.  Same  as  step  1  to  step  3  of  Procedure  A. 

2.  Same  as  step  4  of  Procedure  A. 

3.  Close  valve  VIO,  open  valve  V7  and  Vll  and  change  tempera- 
ture to  the  next  higher  value  for  one  hour  and  repeat  step  2. 
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Figure  B-1.  Pressure  Correction  for  the  Heleod  Gauge  Reading 
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4.  Repeat  step  3  until  460  C  is  reached,  then  cool  to  300°C. 

5.  When  the  temperature  reaches  300  C,  adjust  the  pressure 
to  the  next  higher  value;  also  adjust  bleed  rate  to  the  desired 
value. 

6.  Repeat  step  2  to  step  5  until  the  pressure  of  2310  is 
reached . 
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APPENDIX  C 
FLOW  PATTERN  IN  THE  DIFFUSION  CELL 

The  assumption  that  perfect  mixing  occurs  in  the  diffusion 
cell  will  be  shown  to  be  invalid  for  the  case  when  Equation  (115) 
applied  to  the  experimental  data. 

Assuming  that  perfect  mixing  occurs  in  the  diffusion  cell, 
the  material  balance  for  the  apparatus  shown  schematically  below  will 
be: 


N. 


c> 


Where  F  =  Feed  rate/unit  area 
B  =  Bleed  rate/unit  area 
\  "^  Transport  rate  through  membrane/unit  area 
Xp«5^g  -  Mole  fraction  of  hydrogen  in  F,  B  respectively 
P.Pg  =•  Pressure  of  hydrogen  in  F,  B  respectively 
P_,  =  Total  pressure 
The  material  balance  gives 

F  =•  B  +  N. 


Eliminating  F  from  Equations  (C-1)  and  (C-2)  gives 
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(C-1) 
(C-2) 
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Xp(B  +  N,)   -  BXg  +  N^  (C-3) 

But  for  perfect  mixing  and  at  constant   temperature 

1 


\  -  ^>B-p  (-  i  ) 


i.         1 

=  k  Pg  %  k  Pg  (C-4) 


And  from  the  Ideal  gas  law 


'^B-Zn^)   C^)'"^  ^'-'^ 


where  n„  =  Gram  moles  of  hydrogen  in  the  bleed 

n-  =   Total  number  of  gram  moles  of  gases  in  the  bleed 

Therefore 

^A  "  ^  ^Vb^  "  ^   ''b  (*^"^> 

solving  for  x^  gives 
o 

X3  -  k^l  (0-7) 

Substituting  the  value  of  x_  from  Equation  (C-7)  into 

D 


Equation  (C-3)  and  rearranging 


^  /  B  \  .  ''W 


x„  -  1 


"""  7"  K\  )  ^  IT  "=■" 


or 


'"(^) 


"a"-"!!    rr^zj  ■    <■=-" 


Thus  Equation  (C-9)  is  seen  to  differ  from  Equation  (115). 
This  implies  that  perfect  mixing  does  not  occur  in  the  diffusion  cell, 
and  that  the  actual  flow  situation  lies  somewhere  between  plug  flow 
and  perfect  mixing. 
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